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A study has been made of the Nemst effect and the variation of thermoelec- 
tric power unth magnetic field for single-crystal samples of n-type germanium 
of various orientations and impurity concentrations, at fields up to 18,000 
gauss and temperatures from 275° to 60°K and below. Except at the highest 
temperatures, both effects arise predominantly from that part of the thermo- 
electric power which is due to phonon drag. All observations can be quanti- 
tatively accounted for by theory. They yield information about the dependence 
of the relaxation limes for phonon-phonon scattering on the frequency of 
the phonons, and establish with some certainty the conclusion that four- 
phonon collisions are much less important than thrce-phonon collisions in 
the pertinent range of temperatures and phonon frequencies. Auxiliary in- 
vestigations have shown that the quantization of electron orbits in a magnetic 
field has only a small effect on transport properties when the cyclotron level 
spacing is less than thermal energy. The mean free path of electrons is shown 
to be energy-independent, as acoustic-scattering theory predicts. The absolute 
mobility can be predicted to within 10 per cent or better from data on the 
fractional changes of resistance with stress and magnetic field . 

A more detailed summary of the conclusions and implications of the pres- 
ent work is given in Section IX. 

I. INTRODUCTION 

In a recent paper 1 we have presented measurements of the Nemst 
effect and the change of thermoelectric power in a magnetic field for 
n-type germanium of high purity. These two effects were shown to be 



* This paper, though complete in itself, constitutes part II of a study, part I of 
which was presented in Ref. 1. In addition to describing the phenomena and the 
physical principles underlying them, Ref. 1 develops the background of the present 
study and gives details of the measuring techniques. 
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compounded out of electron-diffusion and phonon-drag contributions, 
the latter being the predominant one in the range near liquid-air tem- 
perature. Now the phonon-drag phenomenon 2 — the pushing of charge 
carriers from hot to cold by the asymmetric phonon distribution which 
a thermal gradient produces — depends on the details of both the elec- 
tron-phonon interaction and the processes which scatter phonons. The 
anisotropics of these interactions and their dependences on wave number 
enter in different ways into the several thermomagnetic quantities which 
one can measure. Thus, an analysis of thermomagnetic data should 
make it possible to sort out the different factors involved in phonon 
drag, and to obtain information not previously available about electron- 
phonon and phonon-phonon interactions. In our first paper we used 
qualitative arguments to draw a number of semiquantitative conclu- 
sions about these interactions. Our object in the present paper is to 
make these conclusions as quantitative as possible, by comparing the 
observations with explicit theoretical formulas. 

1.1 Basic Concepts 

Following our previous practice, 1 we shall formulate the theory in 
terms of the Peltier tensor n of the electrons, which is related to the 
thermoelectric power tensor Q in a magnetic field H by 

Q„(H) = IIft ' ( ~ H) ■ (1) 

Here, n is defined as the tensor relating the energy flux F, relative to 
the Fermi level, to the current density j: 

F = n j. (2) 

In terms of the Q-tensor, the change of thermoelectric power in the a 
direction, due to H, is 

AQ = Q Qa (H) - Q aa (0). (3) 

The Nernst coefficient B(H) is defined — for the symmetrical orienta- 
tions considered here — by 

E„ = -BJ3. X V7\ (4) 

where E N is the open-circuit field transverse to the temperature gradient. 
We shall obtain B from ( 1 ) and the relation 

BH=-Q yx , (5) 

which applies if VT is in the ^-direction and H in the ^-direction. 
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The model which we shall use as a basis for our calculations is the 
"electron-group" model described in our previous paper. 1 In this, the 
possible states of motion of the charge carriers are divided into a num- 
ber of groups, g, and the assumption is made that for each such group 
the distribution function over the group is determined by the contribu- 
tion, j g , which the states of this group make to the current density. 
Carriers of all groups are, of course, subject to the same electric and 
magnetic fields. 

For a multivalley semiconductor like n germanium, a natural choice 
is to let a group g consist of the states in an ellipsoidal shell of the Bril- 
louin zone, of energy range e to e + dt, in a particular valley. Each such 
group g has its characteristic Peltier tensor Tl g ; this must be independ- 
ent of H if, as we are assuming, H affects the j ff 's of the different groups 
but not the distribution function for a given j„ . The total Peltier tensor 
of the medium can be compounded out of those of the different groups. 
Explicitly, 

n(H) = s n • d,(H) • e (H), (6) 

where 6 g is the conductivity tensor of group g and p = (S ff d 9 ) _1 is the 
total resistivity tensor. 

For any one of the ellipsoidal shells just mentioned, the Peltier tensor 
n must have the symmetry of the valley, which for n germanium is 
axial symmetry about a [lll]-type direction. Thus, H g is describable in 
terms of its two principal components, IT||(e), along this direction (the 
high-mass direction), and U x (t), normal to it. Each of these components, 
in turn, is a sum of an electron -diffusion term (subscript e) and a phonon- 
drag term (subscript p). Here, n e || = U ex and is a linear function of e, 
while IIp||(e) =^ n pJ .(e). Our task is to relate these functions of energy 
explicitly to the low- and high-field values of the Nernst coefficient B 
and the change AQ in the thermoelectric power. 

1.2 Qualitative Properties of R p \\ and U p± 

As we formulate the detailed theory in the following sections, it will 
be helpful to keep in mind two qualitative conclusions which can be 
drawn from our earlier study. The first of these concerns the anisotropy 
ratio of the phonon-drag Peltier coefficient. It is that 

54^ »1 but <^^17, (7) 

n pJ .(e) w*j.T|| 

where m*n, m* ± are the principal effective masses and th , t x are the 
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acoustic-scattering relaxation times 3 for current contributions in the 
corresponding directions. A large anisotropy of this sort is to be expected, 
since, if all the phonon modes had the same velocity and relaxation time 
T p h , the ratio Il p \\/Il p± would be the ratio of the rates of crystal-momen- 
tum loss to the lattice for unit currents in the high-mass and low-mass 
directions, and this is just m*i\Ti./m* ± T\\ . 

Our second conclusion had to do with the energy dependence of n P || 
and IIpj. . All the quantities which we have measured correspond to 
various kinds of weighted averages of IT M and II j. . A typical such aver- 
age of IIp|| and IIpj. , considered as a function of energy e, decreases as 
e increases, at a rate distinctly slower than e~\ This was shown to cor- 
respond to a frequency dependence of the relaxation time t p h of the 
phonon modes intermediate between uT 1 and uT 2 and nearer to the 
former. 



1.3 Program 

In Section II and Appendix A we shall discuss critically the adequacy 
of the electron-group model, reporting several auxiliary experimental 
and theoretical investigations which bear on this question. After estab- 
lishing that the model should be quite good for pure n germanium at 
liquid-air temperature and above, we shall describe in Section III the 
procedures used for deriving the formulas for B and AQ. The explicit 
derivations will be given in Appendix B in a form applicable to any 
multivalley medium (e.g., n silicon) ; the formulas we need for the pres- 
ent application will be summarized in Section IV. In Section V we shall 
present the raw experimental data and discuss the corrections which 
need to be applied to them before comparing them with the theory. In 
Sections VI and VII we shall make the comparison for temperatures 
in the range 60° to 131°K, with the object of deducing from it as much 
information as possible about the functions n P ||(e) and n pJ .(e). In 
Section VIII we shall discuss various further observations, such as 
the behavior of B at high fields and the thermomagnetic effects above 
200°K. These observations, which for experimental or theoretical reasons 
are less amenable to precise analysis, still are in good accord with the 
conclusions previously drawn, and allow them to be slightly extended. 
Section IX will summarize the arguments and give a few remarks on 
the significance of the conclusions for the theories of electron-phonon 
and phonon-phonon interactions, topics which we hope to treat more 
fully in later publications. 
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II. CRITIQUE OF THE ELECTRON-GROUP MODEL 

For a multivalley semiconductor the electron-group model, described 
above and in our earlier paper, 1 is based on two assumptions: (a) that 
the distribution function in any particular ellipsoidal energy shell can 
be approximated by a linear function of the crystal momentum over the 
shell and (b) that this distribution function is determined only by the 
external fields, being independent of the distribution functions of other 
shells. We must first examine critically the adequacy of this model; then 
we shall consider further specializations of it. 

The adequacy of the electron-group model depends on the answers to 
the following questions: 

i. Is it legitimate to use a transport equation in crystal-momentum 
space? 

ii. How accurately can one separate the transport problems for differ- 
ent energy shells? 

iii. How well is the distribution function in a single shell approximated 
by a linear function of the crystal momentum? 

As regards i, there are two things which might undermine the validity 
of a transport equation in crystal-momentum space alone — inhomo- 
geneities and orbital quantization. Distortion of the current distribution 
by inhomogeneities in the specimen may make it necessary to consider 
the transport problem in position as well as in velocity space. This 
trivial but annoying difficulty in the interpretation of experiments will 
be discussed at length by one of the authors in another place; it causes 
the transverse magnetoresistance to fail to saturate as H — » <x , but has 
very little effect on thermoelectric and thermomagnetic properties. The 
reason for the difference is simple. When a current flows, portions of the 
specimen with different carrier concentrations try to set up different 
Hall fields, and can only adjust to one another by a distortion of the 
current lines, which becomes large as H —> °o , Thermoelectric and ther- 
momagnetic effects, on the other hand, are measured at zero average 
current, and set up fields which are much less sensitive to inhomogenei- 
ties than the Hall field is. The fields also remain finite as H — * oo , whereas 
the Hall field does not. For this reason, and because measurements on 
different pairs of side-arms showed the gross homogeneity of our speci- 
mens to be quite good, we believe the effect of inhomogeneity on our 
results to be negligible. 

A less trivial objection to the use of a transport equation in crystal- 
momentum space is the fact of orbital quantization, which is known to 
make the high-field magnetoresistance qualitatively different from the 
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predictions of a transport theory formulated in crystal-momentum 
space. 5 " 11 It can be argued theoretically that this is not serious until the 
spacing, fuo c , of the cyclotron levels becomes comparable with kT; how- 
ever, for n germanium with 18,000 gauss along a [100] axis, hw c corre- 
sponds to kT at about 25°K, and one may wonder whether this is negli- 
gible for experiments at 77°. Fortunately, theoretical and empirical 
evidence, which we have gathered together in Appendix A, indicates 
that for our specimens the effect of orbital quantization on the resistance 
and thermoelectric power is probably not more than 5 per cent or so at 
18,000 gauss and 77°, and that for many orientations it is probably con- 
siderably less. The main arguments can be summarized as follows: 

1. The theory of longitudinal magnetoresistance, when interpreted 
with allowance for the energies of the phonons which do the scattering, 8 
predicts only rather small departures from the unquantized theory as 
long as ho) c < kT. Observations 12 support this conclusion. 

2. Both for longitudinal and transverse cases the phonon-drag part 
lip of the Peltier coefficient is proportional to the product of the resistiv- 
ity, the average relaxation time of the phonon modes which scatter the 
electrons, and the average squared velocity of these modes. As the effect 
of orbital quantization on the average relaxation time and velocity is 
probably rather less than on the resistivity, the quantization effect should 
be of similar magnitude for electrical and for thermoelectric measure- 
ments. 

3. Although the theory for transverse magnetoresistance is more com- 
plicated than for longitudinal magnetoresistance, a combination of rea- 
sonable arguments with empirical evidence yields a rough upper bound 
for the orbital quantization effect. Contradictory approaches to the 
theory have been elaborated for the case of acoustic scattering by 
Klinger and Voronyuk, 9 * Argyres 10 and Wolff; 11 the resulting formulas 
can be evaluated for high fields (ftco c ^ kT), but become very cumber- 
some at lower fields. A plausible presumption is that for constant carrier 
concentration and in the limit of small scattering (Hall angle practically 
t/2) the departure from the resistivity given by the older transport 
theory varies at least as fast as 1io) c /kT, as long as this quantity is rather 
less than unity. The size of the factor of proportionality cannot be esti- 
mated reliably from the nonsaturation of transverse magnetoresistance, 

* These authors evaluate their expressions only for the limiting case of fields 
so large that the energy of an average phonon causing transitions of the electrons 
is » kT. However, it is easy to show that, for the case h<a c /kT y> 1 and phonon 
energies negligible, this approach yields a magnetoresistance that goes very nearly 
as H 2 . (Note added in proof: We, are indebted to E. N. Adams and T. Holstein 
for communicating to us an as-yet-unpublished analysis that shows the approach 
of Ref. 9 to be the correct one.) 
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since this can be shown to arise in large part from other effects. However, 
the factor can be estimated from thermoelectric data at liquid-hydrogen 
and intermediate temperatures. 

4. The fair agreement of theory and experiment on high-field trans- 
verse magnetoresistance would be spoiled if orbital quantization in- 
creased the resistivity by any sizable amount. The resulting discrepancy 
would be hard to explain away, as most perturbing effects (inhomoge- 
neities, surface conduction, etc.) also increase the magnetoresistance. 

Turning to question ii at the beginning of this section, we must con- 
sider scattering processes which change the energy. In highly doped ma- 
terial the most important processes of this sort are electron-electron col- 
lisions. There is also another effect : at carrier concentrations high enough 
to make the low-frequency phonons have an effective drift velocity com- 
parable with that of the electrons — the so-called saturation effect 2 on 
lip — the contribution of any electron group to II is reduced by the un- 
balance of the distribution function of the phonons with which it inter- 
acts. Since this unbalance is influenced by the interactions of these 
phonons with all other electron groups, it is no longer possible at such 
concentrations to separate n • j into a sum of n„ • j„ . At low carrier 
concentrations, however, both these effects become negligible, and the 
only type of energy change on scattering which must be considered is 
that due to the finiteness of the energies of the phonons. Only acoustical 
phonons need to be investigated, since scattering by intervalley or non- 
polar optical modes occurs with almost equal probability to all final 
states in an energy shell, and therefore is describable by a relaxation 
time dependent only on the initial state. 

The most detailed study which has appeared on the effect of acoustic- 
phonon energies is that of Dorn. 14 He estimated the effect on the distri- 
bution function for a simple-model semiconductor in an electric field, 
using both a high-temperature expansion method and a variational 
method. He found that, as one expects, the fractional alteration in the 
mobility due to the phonons having finite rather than zero energy is 
very small in the temperature range of interest to us. The alteration is 
of the order of m*e/kT, where m* is the effective mass and c the velocity 
of sound. With a reasonable average m* for n germanium, m*c/k is only 
a small fraction of a degree. However, Dorn's expressions show that the 
fractional alteration in the current of a shell of energy e is of order m*c 2 /e 
for e <$C kT but ~2> m*c. Since for pure lattice scattering the low-field 
magnetoresistance and AQ effects are dominated by the behavior of low- 
energy carriers, the effects of the phonon energies may be expected to 
be considerably more serious for these than for the mobility in the ab- 
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sence of a magnetic field. A quantitative estimate of these effects in a 
magnetic field and with anisotropic masses would be of some interest. 
However, it seems likely that they are still small. The fractional altera- 
tion in the H 2 term in the current of a shell should still be of order 
w*c 2 /e, although perhaps with a slightly larger coefficient than for Dorn's 
case. Since, for a shell c to e + de containing dn carriers, this H 2 term 
goes as er 3 dn « de for e « kT, the phonon-energy correction should 
involve / de/e. Thus, we guess that 

. , m*c , kT , -. 

fractional correction to magnetoresistance tt -y^ In —^- 2 . (o) 

Since the phonon-drag Peltier coefficients of the different shells do not 
vary much with energy, (8) should apply to the thermoelectric power 
also. For n germanium in the range near liquid-air temperature, the 
value of (8) is at most a few per cent. 

We can now give a prescription for answering question iii, the last of 
the three posed at the start of this section. This question had to do with 
the legitimacy of approximating the distribution function in each energy 
shell by a linear function of crystal momentum. Since questions i and ii 
have by now received favorable answers, we need only to solve the trans- 
port equation for an energy shell under the influence of arbitrary electric 
and magnetic fields, and see how well the current and Peltier-heat con- 
tributions from the accurate solution agree with those from the electron- 
group approximation. Now the latter approximation amounts to retain- 
ing only I = 1 spherical harmonics in the distribution function in a new 
set of variables which take the ellipsoidal energy surfaces into spheres 3 
— variables which we have used for other purposes in Appendix A. The 
terms in the scattering operator which mix I = 1 with higher I values — 
call these Su — can be shown 3 to affect the electric current only to the 
second order, i.e., | Su | 2 . They may, however, affect the Peltier flux to 
order S u Ti , where Ti measures the Peltier flux which would be pro- 
duced by a distribution with the I value in question. Now the longitudi- 
nal and transverse branches individually contribute quite anisotropic 
scattering, but their combined effect in germanium is almost isotropic, 
so the Si/ for I t^ 1 are small. But since the relative contributions of the 
two branches to the phonon-drag Peltier effect are not known, it is quite 
conceivable that the S u Ti for I 7* 1 are not negligible, even though the 
| Su | 2 are. A detailed investigation of this question using deformation- 
potential theory has been made for the case I = 3 and, to some extent, 
for larger I, and will be reported elsewhere. 15 The result is that, for ger- 
manium, both the Su and the Ti are quite small, and that the retention 
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of only the / = 1 term in the distribution function should be a very good 
approximation for all conduction and thermoelectric phenomena, with 
or without a magnetic field. 

III. PROCEDURE FOR COMPUTING B AND AQ 

Having satisfied ourselves that the electron-group model described in 
Section I is a good approximation for the cases we wish to analyze, let 
us now use this model to evaluate explicitly the expressions for B and 
AQ which were given there. In both cases the task is to evaluate (6). 
This determines the Xernst coefficient B, or the Nernst field (4), via 
(1) and (5). Similarly, it determines AQ via (1) and (3). Since each 
group g consists of the set of states in the energy range e to € + dt in 
some valley i, it is appropriate to replace the conductivity 6 of this 
group by the infinitesimal dd u> , and to rewrite (6) in the form 

n(H) = 2,T/n <0 (e) .dd (i) (H)] - e (H), (9) 

where the tensor n (,) (e) has principal components ily(e) and n x (e) in 
the principal-axis system of valley i, and where the integration is really 
over energies e. The tensor rfd (,) (H) can be obtained by solving the 
transport equation relating the isothermal current d] ' in the energy 
shell to the electric field E producing it, and expressing the solution in 
the form 

dj U) = d6 (,) ■ E. 

For the case of nondegenerate statistics, which we shall assume 
throughout, the transport equation for dj (,) is easily written down, in 
the electron-group approximation, by equating the rates of gain and 
loss of crystal momentum in the shell. Anisotropic scattering can be 
taken into account 3 by assigning to the shell a relaxation-time tensor t, 
with principal components T||(e), rj.(e) along and normal to the axis of 
the valley. The resulting transport equation — in which, for simplicity, 
we omit the valley suffix i — is 3 






m* • d] ± Mj) dj X H = e- [£) E dn, (10) 



where dn is the number of carriers in the shell, m* is the effective-mass 
tensor and (e) = (3kT)/2 is the mean energy, and where the upper sign 
is for electrons, the lower for holes. For small H this can be solved by 
iteration to give the contributions (48) and (49) of Appendix B to the 



Pa? 



(1) 



PaP 



(+1) 
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conductivity tensor d6. As H — > » the part of E normal to H becomes 
interpretable as the Hall field, and the asymptotic value of the part of 
d] normal to H can be determined by taking the cross product of (10) 
with H and throwing away the first term on the left. This part of dj is 
of order fl -1 ; the H -1 part of dj parallel to H can be determined from 
it by dotting (10) with H. Higher powers in the expansion of dj in H~ l 
can be obtained by further iteration. The resulting contributions to da 
are given in (68) through (72) of Appendix B. 

The resistivity tensor p(H) occurring in (9) is just the reciprocal of 
2,- / dd w (H). A few properties of g and dd (i) are worth listing for refer- 
ence. With superscripts 0, 1, 2, • • • to denote coefficients of 1, H, H , 
in the expansion of p or rfd in powers of H, and +1, », — 1, • • ■ to denote 
coefficients of fl, 1, fl" 1 , • ■ • in the expansion in powers of H~\ we have, 
for a cubic crystal, 

= -R(0)^afiyH y , (11) 

= -fl(«>)2 7 W7 7 , (12) 

where R(0) and R( °° ) are the limiting values of the Hall constant at 
H = and H = °°, respectively, and where 8 a p y = ±1 when afiy is an 
even (odd) permutation of 123, zero otherwise. For n even (odd), p (n) 
and dd {n) are symmetrical (antisymmetrical). The symmetrical tensor 
p (o0) is finite and positive definite, but the only nonvanishing component 
of dd (oo) is the HH component. 

The procedure of Appendix B is thus to get the high- and low-field 
expansions of dd (i) (H) by solving (10); to use these, when necessary, 
to get the corresponding expansions of p(H); to substitute into (9) to 
get the expansions for various components of n(H) ; and thence to com- 
pute the high- and low-field B and AQ for various orientations. The re- 
sults, containing integrals over II||(e) and n 1 (e), can be expressed in a 
great variety of forms. Although in Appendix B expressions are obtained 
for B and AQ valid for aribitrary functional forms of TT 1( , n x , and for 
valleys of either the [111] or [100] types in a cubic crystal, these expres- 
sions are cumbersome in their most general forms. These forms involve 
Maxwellian averages of quantities eIl||, J ./ (n> (T|| , t a ), where / (n) is a homo- 
geneous function of degree n in its arguments. These forms therefore 
simplify greatly if the ITs and t's are each proportional to some power 
of the energy e. Now n e || = lT eJ . is linear in e, while our previous analy- 
sis 1 indicates that n„n and II pJ . probably vary at a rate between con- 
stancy and e _i . Moreover, for ideal acoustic scattering and negligible 
phonon energies tj and t± are a c - *. Thus, specialization of the formulas 
to power-law dependences of the IT's and t's can give a fair picture of 
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the behavior of n c and n p separately, although not as accurate a picture 
as we shall need for our final analysis. We shall give such a specialization 
in the next section, before discussing more complete formulas. 

IV. COMPILATION OF FORMULAS 

Tables I through IV give formulas for high- and low-field B e , B p , 
AQ, and AQ P , as obtained by the procedure just described for various 
cases corresponding to energy-independent anisotropics. The formulas 
are for [111] valleys and are specialized to the cases II P || iX a: e or inde- 
pendent of e, and T||, x cc € ~' or independent of e. Comparison of the latter 
two alternatives makes possible a crude estimate of the effect of impur- 
ity scattering. Although these special cases do not give as accurate a 
representation of the data as we shall ultimately wish to use, they do 
give a fair representation. They are worth listing because their compara- 
tive simplicity makes it easy to see how sensitive the various measurable 
quantities are to the assumed anisotropics and energy dependences. All 
the tables contain references to equations of Appendix B. Most of these 
equations are more general expressions for the quantities tabulated, valid 
even when the anisotropics are not energy-independent. Also, they are 
applicable to multivalley band structures different from that of ger- 
manium. 

The formulas for the low-field Nernst coefficient (see also Equation 
B.4 of Ref. 1) are most conveniently expressed in terms of the dimen- 
sionless coefficients f e and f p , defined by 

WW (13) 



b p = r„ i Q f 



(?)■ 



where, as always, B e and B p represent, respectively, the electronic and 
phonon-drag contributions to B. Table I gives the expressions obtained 
in Appendix B for £% and f p for the case of energy-independent aniso- 
tropics. Note that (1 + f p ) is the product of a factor dependent only 
on the energy variations of the n's and t's by a factor dependent only 
on their anisotropics. The anisotropy factor is unity when the tensor 
m* • -c -1 is isotropic, and likewise when the tensor n p0 is isotropic 
(n p || = n pJ .) ; the energy factor is unity if either -c or U PB is independent 
of energy. Fig. 1 shows the dependence of the anisotropy factor in 
(1 + f p ) on p = n p ||/n px . Curves are drawn for several values of the 
anisotropy w of m* • «c~ ! , chosen to encompass the ranges likely to occur 
for germanium and silicon. 
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Table I 
Values of the dimensionless coefficients f „ , f p , related by (13) to the 
electronic and phonon-drag parts of the low-field Nernst coefficient, 
respectively. The expressions given apply when the ratios w = T\(m ± */ 
t x W||* and p = ILpu/IIpx are independent of energy. They are valid for 
any type of valleys in a cubic crystal (e.g., n silicon or n germanium). 
Angular brackets represent Maxwellian averages, i.e., averages with 
weight e 1/2 exp (-e/kT). For t, one may take either t {1 or t x , and, for 
lip , either n p| | or n pJ . . 



General formula 



Reference, Ap- 
pendix B 

Value for: 
t oc 6 -i/2, 

t independent of 
e, any H„(«) 

Any t (e), ^.in- 
dependent of e 



3 r w _ &v 

<er%> <e)(er>_ 

(57), (50) 



2l 



i + r P 



<dfri*Xer) j~ (2 + u>)(l +w + pw) l 
<*»»<«•*> L (1 + 2u>)(2 + pw) J 

(57), (51) 



4/ir [as above] 
1 [as above] 
1 [as above] 



Table II gives the part of BH going as H' 1 when H — * =o , for certain 
special directions of H and j, again as obtained by specializing the for- 
mulas of Appendix B to the case of [111] valleys and energy-independent 
anisotropics. Note that at high fields the electronic contribution to the 
Nernst field is anisotropic and depends on the anisotropy w of m* • <T , 
whereas at low fields it does not. Fig. 2 shows the dependence of the 
phonon-drag contribution on p = ll pVt /n pJ . , for two assumptions about 
the energy dependence of H ph± . Several features of the graphs and for- 
mulas are worth noting: 

i. The high-field Nernst coefficient is very much larger when H is in 
an [01 l]-type direction than in an [001]-typc. This could have been antici- 
pated from the familiar difference in rapidity of saturation of the Hall 
coefficient in these two directions, and the close relationship of the 
Nernst and Hall effects. 1 The cause in both cases is the fact that, for 
H along [Oil], two of the four valleys have a high-mass direction exactly 
normal to H, and so require large fields to obtain Hall angles near ir/2. 

ii. The theoretical behavior of the Nernst coefficient differs from that 
of the Hall coefficient in that, when H is along [Oil], the value of B is 
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Fig. 1 — Dependence of the anisotropy factor in the expression of Table I for 
1 + f ,, on p = Ilpn/IIpj. , for various values of w = m*xTJm*\\Tx . 

slightly different for VT along [100] and [Oil]. Equality of Hall coeffi- 
cients for these two directions of j is required by the Onsager relations 
(principle of microscopic reversibility), but equality of the Nernst co- 
efficients is not. For the phonon-drag B p the theoretical difference is 
very small when n p || y> II pJ . , but for the electronic B e it is sizable, e.g., 
a factor 1 .62 for acoustic scattering and id = 0.06. 

iii. From Tables I and II it follows that the ratio of B p to B c is not 
usually the same at high fields as at low. For most specimens, therefore, 
there should be a range of temperatures — near that for which B p + 
B c = — where the sign of the Nernst voltage will change with increase 
of H. Substitutions in the formulas give the result that, for any values of 
the anisotropics in the neighborhood of those obtaining for pure germa- 
nium (p zz 10, w « 0.06), the ratio | B p /B e | is greater at high fields than at 
low fields whenever H is along [Oil]. When H is along [001], this ratio 
is less at high fields than at low if U p \, x increase with decreasing energy, 
but becomes the same at high and low fields if IIp^j. are independent of 
energy. 

The electronic part of the thermoelectric power, Q C (H), behaves very 
simply us H —^ co . When H and j are parallel, the energy distribution 
of the current is the same at H = » and H = 0. So, f or H || VT, 



Q.(«) = Q.(0) ==*=£(- 



m 






(14) 
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Fig. 2 — Dependence on p = II W |/n p j. of the phonon-drag part of the Nernst 
coefficient at high magnetic fields, for two orientations of a cubic crystal with 
[111] valleys. All curves have been computed assuming?/; = >n* X T\\/»i* U T± = 0.06, and 
pure acoustic scattering (r<*e~" 2 ). The full curves apply for n ;) a:e _ " 2 , the dashed 
for lip independent of energy. The curves are fairly insensitive to the energy de- 
pendence of t, but depend strongly on w, the lower curves being for small w almost 
proportional to w, the upper to w~ x . For H along [Oil J and V'/' along [100] (not 
shown) the dashed curve would coincide with that for [011][011], while the full 
curve would differ by only a few per cent over most of the range shown. 
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Table III 
Values of the saturation ratio Q P (H -» eo)/Q p (0), for special directions 
of H and vT. Assumptions and notations are the same as for Table I, 
except that the expressions apply only for [111] valleys (n germanium). 
The energy factor (&,){&)/ {*){&&) has the value 8/(3*0 - 0.849 
when rip « e ~ 1/2 , and the value 1 when either n„ or r is independent of 
energy. 



H, vT 



[001], [001] 
[0111, [011] 

[111], [111] 

[011], any J_H 
[011], [100] 
[001], [Oil] 

[oiu, [nil 



Reference, 
Appendix B 



(78) 
(79) 
(80) 
(81) 
(81) 
(81) 
(81) 



Value, General 



(2 + w) (p + 2 ) 
3(2 + pw) 

(2 + w)(l + w + pw) 
(1 + 2m>)(2 + pw) 

(2 + m>)(6 + P + 2pw) 
(7 + 2w)(2 + pw) 

(eHpXcr) r(2 + w)(l + w + pw) l 
(e){eU,T) L (1 + 2u>)(2 + PW) J 

(dlgXg) f 5 + to + P + 2pw ~| 
(eXeHpT) L 3(2 + pw) J 

Wfrr) [ (2 + w) (p + 2) 1 
<€><en» L 3(2 + pw) J 

<ellp><er> f!3 + 5w + 5p + 4pw ~| 



<«)<6n p r) 



9(2 + pu>) 



where | e b - e f | is the distance of the Fermi level from the band edge, 
the upper sign is for electrons, the lower for holes, the angular brackets 
denote Maxwellian averages, and, for the present case of energy-inde- 
pendent anisotropics, t may be taken as either t\ or r x (only the energy 
dependence matters). When H and j are perpendicular, the energy dis- 
tribution of the current is the same as it would be for H = 0, x = con- 
stant. So, for H _L VT 1 , 



Qe( 



[i 



- -< i • 



+ 



66 — €/ 

eT 



a 



(15) 



Table III gives the formulas for Q„( « ), the saturation value of the 
phonon-drag part of the thermoelectric power, again for [111] valleys 
and energy-independent anisotropics. The ratio Q P (*>)/Qp(0) is the 
product of a factor dependent on the energy variations of n pa and m*--c 
by a factor dependent on their anisotropics. Figs. 3 and 4 show how these 
anisotropy factors depend on H p]l /U p± , again for different assumptions 
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Fig. 3 — Dependence on p = n,,| : /n ; ,j. of the expressions of Table III for Q r (*>)/ 
Qi'(U), for H parallel to V'/' in various directions of a cubic crystal with [111] val- 
leys. The fidl curves apply for w = m*iT tt /m*frx = 0.06, the dashed curves for 
w'= 0.08. 



011 the anisotropy of t. For longitudinal effects ( vT || H) the energy 
factor is unity; thus the ratio for these eases gives information on the 
anisotropy of X\,,„ independently of its energy variation. For the trans- 
verse effects the energy factor goes to unity if either n p0 or -c is inde- 
pendent of energy. 
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Fig. 4 — Dependence on p = n^U p ± of the anisotropy factors in the expres- 
sions of Table III for Q P (°°)/Q P (0), for various orientations of a cubic crystal with 
[111] valleys, when H is normal to vT. The full curves apply for w = wi*j.t,|/ 
m*px = 0.06, the dashed curves for w = 0.08. Note two of the factors are the 
same as for Fig. 3. No curves are given for H || [0111, vT || [111], as this case is 
just one-third the way toward [011][100] from [011][01i]. 
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Table IV, finally, gives the expressions for the low-field limits of 
AQ e /H 2 , AQ P /H 2 . These have a less simple behavior than the high-field 
limits. However, it is noteworthy that, when H is in the [001] direction 
and VT in the [100], the expression for AQ P /Q P involves the anisotropy 
of rip only through the same coefficient f p as occurs in the expression 
(13) for B p . The reason for this has been elucidated in our previous 
paper. 

V. DATA AND CORRECTIONS 

5.1 Measurements and Extrapolations 

Measurements of electrical and thermomagnetic properties were made 
on a number of samples of n-type germanium having various orientations 
and dopings, at a number of temperatures from the liquid-hydrogen 
range to room temperature, and for magnetic fields up to 18,000 gauss. 
The samples were of the usual shape with three pairs of side-arms. Heat 
or electric current was fed in at the ends, the potential drop was meas- 
ured between a side-arm of the first pair and one of the last pair, and the 
Hall or Nemst voltage was measured between the two middle side-arms. 
A full description of the experimental procedures has been given in our 
earlier paper. 

Table V gives all the thermomagnetic quantities measured, as func- 
tions of magnetic field, for temperatures above 60°K; measurements were 
also made at a number of fields intermediate between those listed. Such 
data may be useful for comparison with future theoretical expressions 
for the various quantities in the intermediate-field range. In this paper, 
however, we have developed the theory only for the limiting cases of 
high and low fields, and so the most appropriate quantities for compari- 
son with theory are the high- and low-field limits; these are given in 
Table VI. The limits of the various quantities as H — > were obtained, 
since the quantities are even functions of H, by plotting the quantity in 
question against H 2 , fitting the low-field part of the plot to a parabola, 
and noting its intercept. The average size of the random errors involved 
in this procedure can be gauged from the sample plots given in Fig. 5. 
It will be noted that, at low fields, AQ is a little less accurately deter- 
minable than Ap, because of the smaller signals involved; similarly, the 
errors in B are larger than those in R. For the same reason, the random 
errors in the thermomagnetic quantities become larger, percentagewise, 
at higher temperatures. The extrapolations to H = «> Avere made in the 
same way, by plotting the quantities against \/H" and fitting with parab- 
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olas, as shown in Fig. G. The legitimacy of assuming that AQ and Ap are 
represent able by power series in 1/H 2 is borne out by the calculations of 
Appendix B, within the limitations imposed by the neglect of orbital 
quantization, inhomogeneities, etc. However, these latter factors are 
known to prevent Ap from saturating completely as H — > <» , so the ex- 
trapolations may be slightly different from the value which it is appro- 
priate to compare with the present theory. We shall discuss this further 
below. Random errors in the high-field extrapolations are, however, in- 
consequential, as Fig. 6 shows. 

A general discussion of the sources of experimental errors and the 
probable accuracy of the data has been given in our previous paper, 1 and 
will not be repeated here. However, there are some new observations 
which have an important bearing on the trustworthiness of the results, 
and in the following paragraphs we shall discuss these and some theoreti- 
cal adjustments which should be made to the data in order to make them 
comparable with the theory. We have tried to evaluate the more predict- 
able of these corrections for all the table entries that we shall use in the 
analysis in Sections VI and VII, and we have given the sum of these 
corrections in parentheses after the quantity to which they apply. 

5.2 Surface Damage 

In the attempt to track down empirically the magnitude of the effect 
of orbital quantization (an investigation which we have summarized in 
Appendix A) we encountered a stumbling-block in the high value of the 
infinite-field extrapolation of the resistivity normal to // when // is in a 
[100] direction. At 77°K, for example, the originally measured Ap/p for 
Sample 600 extrapolates to a value 0.5 unit above the value predicted 
by the electron-group model. This is many times larger than the ex- 
pected effect of orbital quantization as deduced from other evidence. 
However, such an increase could occur if a partial short-circuiting of I he 
large Hall field were effected by surface or dislocation conduction. A 
surface contribution to the conductance of only a few parts in a thousand 
would suffice. To check this, some of the measurements on this and other 
samples were repeated after etching off about a mil of the original sand- 
blasted surface with superoxol. For the case mentioned (Sample 606 at 
77°), two-thirds of the original excess of observation over theory disap- 
peared; the full excess reappeared on sand-blasting again. Thus, it is 
clear that reliable transverse magnetoresistance results cannot be ex- 
pected with sand-blasted surfaces. 

Theoretically, one would not expect most of the other tabulated quan- 
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Fig. 5 — Sample plots of low-field Nernst, magnetoresistance and AQ data 
against the square of the magnetic field H showing extrapolations to zero field. 
Here AVq is the change in the measured thermoelectric emf in volts between the 
end-arms, due to H, at constant temperature gradient. The field H is measured in 
gauss, and B is in volts/degree gauss. 



tities to be nearly as sensitive to surface effects as is the high-field trans- 
verse magnetoresistance, since only in the latter case is a large Hall field 
present. The high-field Nernst coefficient (see Section VIII) is the one 
exception, since its approach to zero makes a small Hall field due to 
longitudinal counter-currents amount to a large percentage error. The 
measurements confirm the expectation that the effects on the other quan- 
tities should be small. However, though small, the changes found in the 
apparent mobility, etc., were often appreciable. The effects seem to in- 
volve a complex interplay of reduced mobility in the damaged surface 
layer and conduction with little or no Hall effect; moreover, the latter 
conduction seems to be extremely anisotropic, since the large excess of 
magnetoresistance was found only when there were [100] surfaces normal 
to the magnetic field. 

Table VI lists, along with the original values for sand-blasted surfaces, 
the values (identified by asterisks) appropriate to etched surfaces, for 
all cases for which the latter were measured. 
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5.3 Side- Arm Corrections 

The raw data given in Tables V and VI were obtained from the meas- 
urements on the assumption that the flow of heat or electric current 
down the sample was purely one-dimensional and that the side-arms 
acted merely as infinitesimal probes with which potentials could be meas- 
ured without disturbing this flow. This assumption is not quite correct, 
for the width of the side-arms (0.06 cm) was not negligible compared 
with the width (0.15 cm) of the current-carrying portion. Currents 
passed through such a sample must bulge out a little into the side-arms, 
in a way which will be altered by a magnetic field. Moreover, the aniso- 
tropy which a magnetic field introduces into the thermoelectric power 
can cause circulating currents to flow near the side-arms in the presence 
of an electric field, even though no net current flows along the specimen. 
These effects can become very serious when a large magnetic field is ap- 
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Fig. 6 — Sample plots of high-field magnetoresistance and AQ data against the 
inverse square of the magnetic field H, showing extrapolations to infinite field. 
Notation and units are the same as for Fig. 5. 
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plied in the direction of the side-arms. For example, in indium antimo- 
nide magnetoresistances of several times the ideal value have been 
observed for this geometry. 16 For this reason, all the transverse measure- 
ments reported here were made with the magnetic field normal to the 
plane of the arms, and we shall show that, for this case and for all longi- 
tudinal and zero-field cases, the side-arm errors are small and can be 
approximately corrected for. 

A number of statements about the effects of the side-arms can be made 
from simple theoretical considerations. For example, it is obvious that 
the current distortions have no effect on the measured values of the 
zero-field thermoelectric power or of the Hall coefficient. The conduct- 
ance of the sample is, of course, higher with the side-arms than without ; 
measurements on scaled cutouts of stainless-steel sheet showed the dif- 
ference of true and apparent resistivities in the absence of a magnetic 
field to be about 2.G per cent for our geometry. However, it can be shown 
fairly easily that the measured transverse magnetoresistance Ap/p, and 
likewise the transverse AQ, are unaffected by the presence of the side- 
arms when H is in a [100] direction ; this in variance arises from the two- 
dimensional nature of the potential distribution and the isotropy of p 
and Q in the plane normal to H. When H is in a [110] direction the aniso- 
tropics of Ap and AQ cause a difference between the apparent and true 
Ap/p and AQ, which for transverse cases can be expressed in terms of the 
anisotropics and the correction to the zero-field resistivity. When H is 
longitudinal, more complicated effects arise, due to the Hall fields where 
the current lines depart from longitudinality at the bases of the side- 
arms. However, for almost all the cases listed in Table VI there exist 
scaling relations which allow the side-arm corrections to be expressed, 
at least approximately, in terms of the correction to the apparent con- 
ductivity at // = 0. 

Several experiments were performed to verify that these effects were 
small and to evaluate the corrections for the two cases (high-field trans- 
verse AQ with H in a [110] direction) for which the scaling relations do 
not suffice for a theoretical prediction. These experiments consisted in 
removing the middle pair of side-arms of some of the specimens and no- 
ting the changes in the apparent values of Ap/p and AQ. The changes due 
to removal of all side-arms would, of course, be twice as great. Table VII 
summarizes the theoretical and experimental results for the corrections 
to be added to the entires of Table VI. In no case are these greater than 
about 4 per cent. 

5.4 Gross Inhomogeneities 

Gradients of impurity concentrations can greatly falsify the results. 
For example, the apparent Hall constant depends on the carrier concen- 
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Table VII 

Approximate correction factors by which the entries of Table VI should 
be multiplied to give the "true" quantities as they would be measured 
on a specimen without side-arms. Except for the AQ'a, the factors are 
independent of temperature, or nearly so. In the last two rows, the 
factors given apply in the neighborhood of liquid-air temperature. 



Direction of j or VT 
Direction of H 

Quantity: 



ft (£) 



lim 



>(//) 



h— p(0) 
B(0) 

lim jvr, 

lim AQ 



[100] 
[001] 


[100] 
[0111 


[H0| 
[110] 


[110] 
[110] 


[100] 
[100] 


1.026 


1.026 


1.026 


— 


— 


1.000 


0.987 


1.006 


0.965 


0.996 


1.000 


0.994 


1.008 


0.974 


0.974 


0.974 


0.974 


0.974 


— 


— 


1.000 


0.975 


1.009 


1.014 


1.013 


» 1.000 


« 1.000 


« 1.000 


1.000 


1.000 



[111] 
[111] 



0.974 



1.000 



tration between the middle side-arms, while the apparent conductivity 
depends on the carrier concentration over the whole length; if the con- 
centration varies, the apparent Hall mobility will be incorrect. Most of 
the samples were checked for homogeneity by measuring Hall constant 
on the end arms as well as on the middle pair. All three measurements 
usually agreed to within 1 per cent. Inhomogeneities can also cause meas- 
urements designed to give quantities even in the magnetic field (Ap, AQ) 
to contain contributions from effects odd in the field (R, B), and vice 
versa. This effect occasionally showed up as an asymmetry in the meas- 
ured values when the direction of H was reversed. This asymmetry was 
small for all cases listed in the tables and was consistent with the present 
interpretation, e.g., the magnet oresistancc asymmetry could be calcu- 
lated from the Hall inhomogeneity. The entries represent averages of 
measurements for H and — H. 

5.5 Value of the Temperature Gradient 



The measurements tabulated in Tables V and VI were all taken with 
the special-purpose apparatus described in our earlier paper. 1 With this 
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apparatus the temperature gradient used for the thermomagnetic effects 
was not measured directly, but was computed from the thermoelectric 
voltage developed between the potential leads at H = 0, together with 
the thermoelectric power Q of the specimen. The latter was computed on 
the assumption that, for a given temperature, the only source of varia- 
tion of Q from specimen to specimen is the variation in the height of the 
Fermi level associated with the different carrier concentrations. In other 
words, the Q for any sample was computed from values previously meas- 
ured on other samples by adding a term equal to 86 /zv/degree times the 
logarithm of the ratio of carrier concentrations. This will not be correct 
if impurity scattering affects Q appreciably. 

Previous empirical evidence suggested that, for fairly pure material, 
Q is almost unaffected by impurity scattering. The same conclusion fol- 
lows from a theoretical analysis which we shall carry out in Section VII. 
When T ^ 77°K, it turns out that no correction is needed to the data 
for Samples 603, 604, 606 and 610, on which most of the analysis of the 
next two sections will be based, nor for Sample 595. The corrections 
should be small for these samples at T £b 60°K, and for Samples 576A 
and 580 at T ^ 77°. For Sample 601 the correction begins to be percep- 
tible at 94°, and for Sample 596 at 87° it amounts to about 8 per cent of 
the quantities tabulated. 

5.6 Boundary Scattering 

The final — and largest — correction we shall consider has to do with 
the scattering of phonons from the boundary of the specimen. We wish 
to interpret the data in terms of a theory of conduction in an infinite 
homogeneous medium. The side-arm correction has eliminated effects 
of the shape of the sample, but we have yet to correct for the finiteness 
of its dimensions. 

We may start from the observation that phonon-drag effects are due 
to the anisotropy of the low-frequency part of the phonon distribution. 
This anisotropy, due to the temperature gradient, has a certain value at 
depths below the surface that are large compared with the mean free path 
of a phonon. However, immediately below the surface the anisotropy is 
only about half as great, since that half of the phonons that comes from 
the surface has no tangential anisotropy if the surface scatters diffusely. 
Thus, one may say, approximately, that the tensor Q p has one value at 
depths greater than some value A, of the order of a phonon mean free 
path, and has half this value at depths less than A. The effect of this al- 
teration is easily seen to be to make the effective Q P (H) less than the in- 
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terior value by the ratio 

Qp =l- 2 -£, (16) 



Q p (ideal) A» 

independently of H, where A is the cross-sectional area of the sample. 
The effect on the low-field Nernst coefficient can be calculated similarly; 
it is 

B " -1-T. (17) 



B p (ideal) t 

where t is the thickness of the sample normal to its plane. We shall take 
A as roughly independent of the orientation of the surface and equal, in 
millimeters, to 0.037 (77° /T) ' , the best value currently available from 
experiments on the variation of Q p with diameter near liquid-air tem- 
perature. 2 

In parentheses after each entry in Table VI we have given the sum of 
the side-arm correction calculated from Table VII and (for the thermo- 
electric quantities) the temperature-gradient correction and the bound- 
ary-scattering correction computed from (16) or (17). This total correc- 
tion is to be added to the table entry to get the ideal value appropriate 
for the comparison with theory. 

VI. ANALYSIS 

At each temperature and for each impurity content we have 10 or 11 
measured quantities which depend on the partial Peltier coefficients 
II||(e), U±(e) of the energy shells, and among which no relations deriva- 
ble from phenomenological theory exist. These are: Q and B at H = 0, 
AQ(H = oo ) for [100], [110] and [111] longitudinal orientations, one 
transverse AQ(H = « ) for H || [100] and two such for H || [110], and 
two (or three) constants describing the low-field AQ behavior. (There 
are three phenomenologically independent constants of the latter type, 
but there is one relation between them which follows merely from the 
band structure and the electron-group approximation. 1 ) We shall now 
derive what information we can about n p ||(e) and H p± (e) by fitting the 
theory of Section IV and Appendix B to the data. 

We shall first note briefly that all these data can be approximately 
fitted by the specialized formulas of Tables I, III and IV, with three 
adjustable constants A» , A ± and n, where 



n„n . = A 



|A (ey 



(592 
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Table VIII 

Comparison of thermomagnetic data observed at 94 °K with values com- 
puted from the formulas of Tables I and III and from the mean moment 
curves of Fig. 7. For the former calculation w = 0.061 and the ap- 
proximately optimum values p = U P JU P± = 9.5, Q p = 697, U p « e -0 - 25 
were used. The corrections in parentheses in Table VI have been included 
in the observed values, and the data for Samples 603 and 610 have been 
corrected to 94°. The data for Sample 606 are starred to indicate that 
this sample had been etched. The units of theAQ's are /iv/degree, those 
of B(0) are 10 _9 v/gauss degree. Acoustic scattering has been assumed; 
correction for impurity scattering would lower the computed B(0) 
slightly. In the first row, Q(0) — Q* is T' 1 times the Peltier heat relative 
to the band edge, i.e., the sum Q p + 172/*v/degree. 







Directions 


Value Com- 




Value Com- 
puted from 
Fig. 7 


Quantity 


Sample 
No. 






puted from 
Tables I 


Observed 
Value 










// 


vr 


and III 




Q* - Q(0) 








869 


869 


863 


-AQ(») 


606 
604 


[001] 


[001] 


1440 


1484* 
1460 


1459 




610 


[HI] 


[111] 


605 


621 


615 




603 


[011] 


[011] 


118 


119 


124 




606 


[001] 


[100] 


91 


144* 


93 




603 


[011] 


[Oil] 


1325 


1191 


1187 




604 


[011] 


[100] 


655 


613 


610 


fi(0) 


606 
603 
604 






22 


31* 

32 

36 


34 



and with the value of w = m* ± T B /m*|rj. given by, say, the high-field 
magnetoresistance measurements. Table VIII gives, in its fifth and 
sixth columns, a sample comparison of observed values with the pre- 
dictions of this simplified theory; low-field AQ values have been omitted 
from the table because, as we shall see in detail in Section VII, they need 
to be corrected for impurity scattering. The extent of the agreement is a 
fairly good measure of the adequacy of the assumed behavior of U p , 
since for most of the quantities listed the electron-diffusion contribution 
is small enough so that uncertainties in its value are unimportant. The 
observations are fitted well enough to engender some confidence in the 
theory, but there are systematic discrepancies which suggest that U v \\ 
and U p± may have different energy dependences that a less restrictive 
analysis could reveal. 
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G.i Analysis in Terms of Moments 

As we have mentioned in Section III, the most general formulas of the 
electron-group theory, as derived in Appendix B, give Q or AQ and B as 
linear combinations of Maxwellian averages of quantities of the form 
cTTjj.a./ '* ( 7-|| , tj.), where f n) is homogeneous of degree n. The coefficients 
contain in their denominators linear combinations of similar Maxwel- 
lian averages, without the n's. Now when n is or 1, corresponding re- 
spectively to the high-field transverse Q and to Q(0), or the high-field 
longitudinal Q, it should be quite a good approximation to take t\\ and 
tj. both proportional to e -2 in the purest of our samples (ideal acoustic 
scattering). The value of w = m* ± T\\/m*\\T± appropriate to these cases 
can be determined from the high-field longitudinal magnetoresistance. 
However, for n = 2 and especially for n = 3 (low-field effects), the 
Maxwellian averages in question become very sensitive to small amounts 
of impurity scattering, and one might question the legitimacy of taking 
the ratio t\\/t ± independent of energy. But since the term containing 
T\\" l T ± "~'" also contains m*\\~'"m* ± '"~", and since II||, X vary only slowly 
with energy, it will always be a good approximation to assume that, 
with angular brackets as usual denoting Maxwellian averages, 

/en -te) m fer>-"(^te))- (i8) 



vte)"'te)"">»»'"('fe)> (10) 

where, for n = 1, 2, 3, w n is determined from purely electrical data. 
Specifically, we can define w„ so that (19) is exact for m = 1; (18) will 
also then be nearly exact for m = 1. For m > 1, (19) and (18) may be 
less accurate, but this inaccuracy will have very little effect on any elec- 
trical or thermomagnetic quantity, because the terms with m > 1, be- 
ing of order w n '", will be very small. 

We shall use the approximations (18) and (19) to express the various 
thermomagnetic quantities in terms of "moments" iV"' and IIj. (n) , where, 
forn = 0, 1, 2, 3, 

n *' x = ~w^r- (0) 

From the formulas of Appendix B we obtain the expressions listed in 
Table IX. Since the low-field Nernst coefficient involves H\\ i± J and 
II||, A (2) , we have listed as the second entry in the first column that com- 
bination of B and Q which involves only the ri|i,j. (2) . Similarly, we have 
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Table IX 
Thermomagnetic quantities in terms of the moments LT||, ± n denned by 
(20), as obtained from the formulas of Appendix B using the approxima- 
tions (18) and (19) and the assumption of n-type material and [111] 
valleys. 



Quantity 



TQ(oo) 



TQ(0) 
r B(0) 



T I -lim 
Mi 



T L(w//c) + 

r-ih 

L H- 

AQ + (BuhW/c) 

(Ap/p) + (w/ff/c) 2 

T r_ nm *Q- 4 

L fl-0 {Ap/p) 



•] 



Directions of 
H and VT 



[Oil] [100] 
[Oil] [Oil] 
[001] any _|_H 
[001] [001] 
[111] [llU 
[011] [011] 



[001] [100] 
[001] [001] 



Expression 



1 ' + 2w nr | Mot^ nV 1 " 



32 + w 



3 (2 + w) 



\ n„(» + I n x <»> 
» n„(o) + 14^.n i <«) 



l+2w 



1 +2w 






1 +2u> 

to 

2 + w 

1 + 2w« 



1 + 2u> 



n„<i> + 



2 + w 



n±m 



n ' w +rr£ ni< " 



" : nv 3 > + — — nx< 3 > 



2 + w 3 



Wu 



1 —IDs 



n.» + 



2 + w a 

1 
1 — w 3 



HiO 



Reference, 
Appendix B 



(81) 
(81) 
(81) 
(78) 
(80) 
(79) 
(50) 
(57) 

(58), (65) 

(58), mod- 
ified (65) 



listed in the last two rows those combinations of the low-field AQ, B, and 
Q which involve only the IT M . x (3) . These equations can in most cases be 
inverted to give the II||, j. (n) in terms of the observed quantities. Thus, the 
first two rows just suffice to determine TT H . j. (0) , the next four rows over- 
determine II||, j. (1) , the next row gives one relation between II|| (2) and U ± l2) , 
and the last two rows just determine TT|j. j. (3) . 

Our procedure will now be to determine the moments U^ <± M from these 
equations and the observational data, and then to discuss the separation 
of each such empirically determined I1||, J . ( " ) into electron-diffusion and 
phonon-drag contributions and to interpret the variation of the latter 
part with n in terms of the energy dependence of n p ||, x . 
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Table X 
Values of w x , w 2 and w 3 used in the moment analysis. The values of Wi 
were obtained independently for each case by fitting the electron-group 
theory to the observed high-field limit of the longitudinal magnetoresist- 
ance. The side-arm corrections of Table VI have been included and, 
whenever possible, data obtained after etching have been used. The 
latter cases are identified by asterisks. The values of w 3 were estimated 
from these W\8 (assumed for the first four samples equal to 0.061) and 
the approximate dependence of w 3 on impurity concentration found by 
Goldberg. 17 The values of w 2 were interpolated between Wi and w 3 . 



Sample 


T, °K 


Empirical wi 


Assumed vit 


Assumed ws 


606 


60° 
63° 


0.061* 


0.0625 


0.064 




77.4° 


0.0615* 


0.0615 


0.062 




94° 


0.061* 


0.061 


0.061 


604 


63° 


0.062 








77.4° 


0.061* 


0.0625 


0.064 




92° 


0.0625 


0.0615 


0.062 


610 


77.4° 


0.061* 


0.0625 


0.064 




92° 


0.061 


0.0615 


0.002 


603 


77.4° 




0.0625 


0.064 




94° 




0.0615 


0.062 


Value adopted for all 


0.0610 






above cases 








601 


94° 


0.069 


0.076 


0.084 


596 


87° 


0.081 






606 


131° 


0.067* 
(0.0595 adopted, 
see text) 


0.0595 


0.0595 



6.2 Choice of the w n 

To begin with, we need values of w, w 2 and w 3 . Magnetoresistance evi- 
dence on w (=Wi) for the present samples is summarized in the column 
labeled "Empirical Wi" of Table X. The good agreement of the Wi val- 
ues from Samples 606 and 604 with that from Sample 610 speaks well 
for the validity of the electron-group model. There seems to be no sig- 
nificant change of W\ with temperature in the range 60° to 94°K. No 
reliable value of Wi can be obtained from Sample 603, since in this orien- 
tation a change of w x from 0.04 to 0.08 only changes the theoretical 
high-field resistance by 0.05 of the zero-field value. The data for 603 do, 
however, agree fairly well with the theory, in that the Ap/p predicted 
from the Wi of the other specimens is within 0.02 or 0.03 of the value in 
Table VI. 
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The empirical w x value for Sample 606 at 131°, though listed, is not to 
be trusted, as a field of 18,000 gauss does not give a very good approach 
to saturation. We have preferred to rely on the value 0.061 found at the 
lower temperatures, decreasing it to 0.0595 to allow for the gradual de- 
crease of w with temperature found by Goldberg, 17,18 a decrease which 
also shows up, though less clearly, in our own low-field data at higher 
temperatures. 

The estimation of w 3 is somewhat less satisfactory. In principle w 3 can 
be estimated from the longitudinal and transverse low-field magneto- 
resistances, since, as we noted at the start of this section, two independent 
phenomenological constants should suffice for the complete description 
of the low-field magnetoresistance in reasonably pure material. However, 
the w 3 values computed in this way for Samples 603, 604 and 610 turn 
out to be quite different — ranging, for example, from 0.056 to 0.073 at 
77°K, all for etched samples — whereas the similarity in impurity con- 
tent between these three samples would lead one to expect similar w 3 's. 
This sort of discrepancy might have been anticipated from the fact 
(Ref. 1, Fig. 13) that the four data for 603 and 604 depart perceptibly 
from the predictions of the phenomenological theory with two adjusta- 
ble constants. The trouble can arise from any of several causes. To com- 
pute w 3 from the data, one needs the longitudinal and transverse mag- 
netoconductances, and the second of these equals the corresponding 
magnetoresistance plus (n H /c) 2 . Now the empirical value of n H , unlike 
magnetoresistance and AQ values, depends on the measured thickness of 
the sample and is affected by errors in this measurement. Also, it can be 
falsified if the carrier concentration between the Hall arms differs slightly 
from the average over the length of the sample. Since a change of mh by 
1 per cent affects w 3 by 0.003, this source of scatter from specimen to 
specimen could well be serious. Moreover, our random errors are large 
enough to affect w 3 perceptibly, and comparison of etched and unetched 
samples has shown that the surface damage effect can alter the empirical 
w 3 by 0.007 or so. 

Fortunately, a study of the fomulas of Table IX, which we shall de- 
scribe below, shows that, once w 3 is specified, the II|| r j. (3) computed from 
the measured quantities is not nearly so sensitive to small errors in the 
latter as we have just found the empirical w 3 to be. This suggests that 
one can get a fairly reliable analysis by disregarding the apparent w 3 s of 
the specimens and estimating the true w 3 for each temperature and im- 
purity concentration in some other way. Now for pure acoustic scat- 
tering w 3 should be practically the same as w x , for which we have ob- 
tained the presumably trustworthy value 0.061 (Table X). Several 
studies of the effect of impurity scattering on magnetoresistance have 
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been published, 1 '' 18,19 and, while the absolute values of w 3 obtained in 
these studies may possibly have been affected by various kinds of sys- 
tematic errors, their picture of the variation with impurity content and 
temperature is undoubtedly correct. From the work of Goldberg and 
of Goldberg and Howard, 18 it appears that, at 77°K, a sample with im- 
purity concentration 7 X 10 13 cm -3 has a iv 3 differing from that of an 
ideally pine sample by only about 10 per cent, or at the very most 20 
per cent.f Thus it is a reasonable guess to say that the corresponding de- 
viation at a concentration of 2 X 10 13 cm -3 is perhaps 5 per cent, i.e., 
0.003 in w 3 . While the error in this estimate may be as large as the esti- 
mate itself, it is surely much smaller than the scatter mentioned in the 
preceding paragraph, and we shall see that this accuracy suffices for our 
present purpose. The last column of Table X gives the values we shall 
adopt for our analysis, obtained by adding an estimate of this sort to our 
«'i . The next to the last column gives the corresponding w 2 's, obtained 
by interpolating between w x and w 3 . 

0.3 Results for II|| (b) and Jl x M 

It is now a straightforward matter to insert the w H of Table X into 
the formulas of Table IX, equate to the empirical quantities of Table 
VI (with the corrections included), and solve for the moments II||, X " 
defined by (20). To make a comparison of the values for different speci- 
mens more meaningful, it is best to compare not the n^j. 00 , which de- 
pend on the location of the Fermi level through n c , but the quantities 
II* - lli,,^" 1 , where 

n* = e S^JH (21) 

c 

represents the part of n„ due to the difference between the Fermi level 
e F and the band edge e fc . The quantities II* — II||, X (,,) thus are positive 
for n-typc material and measure moments of the Peltier heat relative to 
the band edge. 

Fig. 7 shows the results for II* — II ± u,) and II* — Il|| <n) , respectively, 
for a number of the purest samples at several temperatures from 60° to 
131°K. Here the moments for n — 1 were calculated from Q and the 
high-field limit of AQ for H longitudinal. Those for n = '.] were calculated 
from the longitudinal and transverse values of 

,. AQ 

lim — , 

lf-0 II 



t We are indebted to C. Goldberg for having informed us of some of these re- 
sults prior to publication. 
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Fig. 7 — Empirical values of the moments nj. (n) (left) and lV n) (right) de- 
fined by (20), as obtained from the method of analysis described in the text. Note 
that the two scales differ by a factor of ten. The dashed curves at the bottom of 
the chart for nj. (n) show the contributions which IT, would make to n^"' or nj. (n) 
at the various temperatures, if the scattering were purely acoustic. The full curves 
are reference means which we shall use in Figs. 11 and 12. 
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using, if necessary, the pheiiomenological identities mentioned at the 
bottom of Table IV to relate the AQ values for the directions going with 
a particular sample to those for the directions listed in Table IX. For 
n = 2, only one moment can be determined; since II j. " is very repro- 
ducible and varies very little with n, we have chosen to determine n x (2) 
by linear interpolation between TTj.' 1 and IIj. (3) , and then to determine 
II|| (2) for each specimen by the equation specified in Table IX. For n = 0, 
a similar procedure was followed, using a linearly extrapolated value of 
n A ( ; the alternative procedure of combining transverse AQ data on two 
different samples would lead to larger random errors. 

For the 131° data the method just described for determining the mo- 
ments for n = 1 starts to break down, since, as we noted above, the 
high-field limit is not reliable. However, it is probably safe to say that 
AQ( co ) lies in between the extrapolated value of Table VI and a value 
increased from this in the ratio of the Ap( oo ) calculated for the assumed 
Wi to the Ap( oo ) of the table. The first moments obtained from both 
these limits are plotted in Fig. 7 and connected by vertical lines to em- 
phasize the uncertainty. A corresponding though smaller ambiguity re- 
sults for the second moment (II* — II|| <2) ). The second and probably 
more nearly correct assumption gives the higher (II* — H|| ) and the 

lower (n* - n x (1) ). 

Fig. 7 also shows, for comparison, the electronic contribution 

n* - n, (n) , 

as given by an equation analogous to (20) on the assumption t ± <* e , 
i.e., without correction for impurity scattering. The latter correction 
should be very small for these specimens; some estimates of it will be 
given in Section VII. It will be seen that, with our neglect of impurity 
scattering, the electronic contribution is linear in n, and that it accounts 
for the major part of the variation of IIj. (n) with n ; this makes our use of 
linear interpolations and extrapolations seem reasonable. Subtraction of 
the electronic and total moments shows that | U p± . | is almost independent 
of energy, decreasing slightly with increasing energy, and that | n p || | has 
a more rapid, though still modest, decrease. This is a refinement of the 
conclusion of our first paper to the effect that some average of | n p || | and 
| IIpj. | decreases with increasing energy. Since an accurate quantitative 
estimate of the energy dependences of n p || and n px requires a correction 
for impurity scattering, we shall postpone this topic to Section VII. 

The last two columns of Table VIII give a sample comparison of the 
measured thermomagnetic quantities with values predicted from the 
mean moment curves of Fig. 7. The fit is, indeed, significantly better 
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than was obtained in the middle column on the assumption of energy- 
independent anisotropics. 

G.4 Accuracy of the Results 

A few words are in order regarding the sensitivity of these results to 
random or possible systematic errors in the inputs. As one might infer 
from the consistency of the values of uj u and U x l3) for different samples 
and from their regular trend with temperature, these quantities are 
rather insensitive. For germanium with its small w, Table IX shows that 
the quantity n x (l) equals TQ(0) minus a small fraction of the longitu- 
dinal high-field limit of AQ; since the latter term is only about one-fifth 
of TQ p (0), errors in it have little effect. This term is almost proportional 
to w. Thus, for example, it turns out that, at 77°K, an alteration of even 
as much as 10 per cent in the assumed xo, or in the high-field AQ, would 
affect the plotted quantity n* - n x (1) by only 2 per cent of its value. 
The makeup of IIi (8) is similar. It equals TQ(0) plus a number of terms 
involving the Hall mobility and the low-field Ap's and AQ's; the latter 
terms are again only a small fraction of Q„(0), and they depend hardly 
at all on w s . For example, for Sample GOG at 77°K, a change of the as- 
sumed w-i by 10 per cent would alter II* - n x <3> by only one-third of 
1 per cent; changes of 10 per cent in any one of the quantities B, m« 
or the transverse Ap/p// 2 or AQ/// 2 would affect II* - n x (8) by 1 per 
cent or less; a 10 per cent change in the longitudinal Ap/p// 2 or AQ/II 2 
would have a 3 per cent effect. The sensitivity for the other cases is sim- 
ilar. Thus, the values of n x (3> should be quite reliable. 

Turning to the II|| (n) , we find these quantities to be much more sensi- 
tive. The values of n/ , which depend primarily on the longitudinal 
high-field limit of AQ, should be more reliable than those for other values 
of n. Thus, for example, for .Sample GOG at 77°K an error of 5 per cent 
in the high-field AQ would affect II* - II|| 0) by a little over 4 per cent. 
A 10 per cent change in the assumed w would affect this quantity by 
only about one-half of 1 per cent. On the other hand, n/ 3 ' is much 
more sensitive to iv 3 , being in fact very nearly proportional to w 3 ~ . It 
is also fairly sensitive to the [100] longitudinal AQ/// 2 and Ap/p// 2 . For 
example, at 77°K a 10 per cent change in either of the latter quantities 
would change IT* - rT lt <3> by about 9 per cent. However, the sensitivity 
to the other inputs mentioned in the preceding paragraph is small; a 10 
per cent change in any one of them at 77°K would affect n* - II|| <3) by 
no more than 2 per cent. The quantity IT* — II/ 2 * behaves very similarly. 
It is very nearly proportional to wf l ; however, the uncertainties in w 2 
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should be rather less than those in iv 3 . Changes of 10 per cent in B or 
Hh would affect II* — IT|| (2) at 77° by about 5 per cent, and a change of 
3 per cent in the interpolated n x would affect it by 6 per cent. The 
values of II* — II|| (0) , finally, are insensitive to w when derived from 
measurements with H in a [1 10] direction. The sensitivity to the assumed 
high-field transverse &Q is moderate; for Sample 603 at 77°, for example, 
a 10 per cent change in the latter quantity would change II* — II|| (0) by 
perhaps 9 per cent. 

6.5 Oilier Orientations 

Several conceivable calculations have been omitted from the list used 
in the construction of Fig. 7; these correspond to cases where the ob- 
served quantity is very insensitive to the moment which one might wish 
to calculate from it. Thus, as Fig. 3 shows, the longitudinal AQ in the 
[011] direction is very small and insensitive to p (i.e., to n/ ). Similarly, 
Fig. 4 shows that the transverse AQ with H in an [001] direction is small 
and very insensitive to p (or n>| <0) ). Although no reliable points for Fig. 
7 can be computed from these data, the entries in the last two columns 
of Table VIII show that the AQ's for these orientations are in fair agree- 
ment with the predictions of the theory. (It is more significant to com- 
pare values of Q( oo ) = Q(0) -f- AQ than values of AQ itself, since the 
latter quantity could conceivably be of cither sign.) The small discrep- 
ancies remaining arc of the order of magnitude of the expected effect of 
orbital quantization. 

VII. DEPENDENCE OF ITpji , j. OX ENERGY 

The obvious next step is to convert the empirical moments of II* — 
rin.j. plotted in Fig. 7 into moments of — IIp||.j. by subtracting the the- 
oretical moments of II* — II, . One may then hope to determine the 
rate of variation of ^l P \\,± with energy e by comparing the way in which 
these moments vary with n with the /< -variation to be expected for, say, 
n p ||, x « e\ For a specimen completely free of impurity scattering, all 
this could be done very simply indeed. However, as we shall see, the sec- 
ond and third moments, i.e., those with weight factors r x 2 and t ± 3 in 
(20), are extremely sensitive to impurity scattering, since this modifies 
Tx greatly in the low-energy region, where it is large, and since n P ||. x 
also tend to be largest at low energies. Thus, in spite of our efforts to 
avoid the morass of impurity-scattering theory by using very pure spec- 
imens, we must still resort to it if we are to make a complete analysis of 
our data. However, the smallness of the impurity scattering in the sam- 
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pies used for Fig. 7 will turn out to simplify greatly the task of correcting 
for it. We shall therefore digress for a moment to discuss the nature of 
impurity scattering and to establish the validity of the procedure we 
wish to use to take account of it. 



7.1 Impurity Scattering 

The Conwell-Weisskopf theory of ionized-impurity scattering, 20 com- 
bined with various modifications, 21,22 ' 23,24 has had some success in cor- 
relating mobility data and other properties of semiconductors. 25,2 * How- 
ever, a theory of this form can be rigorously justified only in the limit 
of carrier and ion densities rather lower than those normally encountered 
in semiconductors; even for this case, a rigorous treatment requires in- 
clusion of the effects of electron-electron collisions, 27 a refinement which 
has usually been ignored in the semiconductor literature. Under condi- 
tions such as those of our experiments, there are several aspects of the 
theory which become rather dubious, e.g., use of the Born approximation 
in scattering, 23,24 the justification for treating the screened coulomb po- 
tential of an ion as a static scattering potential, etc. Thus, about all one 
can safely say about ionized-impurity scattering in this range is that it 
is something which modifies the distribution function more and more as 
the energy becomes lower. Fortunately, this fact seems to be about all 
one needs to know in order to correlate the various effects of impurity 
scattering with one another, in the range where impurity scattering is 
slight. We shall now try to demonstrate this empirically. 

Although the fact of electron-electron scattering prevents a relaxation- 
time model from being strictly valid for impurity scattering, we shall 
follow custom and assume that, in the range of interest, the effects of 
impurity scattering can be adequately allowed for by modifying the form 
of the functions Ti|(e), T ± (e). The usual assumption takes the form 



= ae' 



+ &e"», (22) 



where the first term represents acoustic scattering, and the second impu- 
rity scattering. This leads to difficult integrals in the various transport ex- 
pressions. Since no real justification can be given for the exact form of 
(22) — except at densities far below ours — we shall try using instead 
another expression, which has the same property of reducing t at low 



* For a general review with further references, see Ref . 26. 
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energies but leads to more tractable integrals, namely, 



r cc e 



D— "(l?)]- (23) 



Here is a parameter which measures the amount of impurity scatter- 
ing. 

7.2 Legitimacy of {23) 

Our argument that (23) is justified will consist of two parts. The first 
will be to show that (23) gives practically the same results as (22) for 
a variety of problems, i.e., that when impurity scattering is just begin- 
ning to be appreciable the exact form of the impurity-scattering law is 
of less consequence that the mere fact that the t's of low-energy elec- 
trons are greatly reduced. The second argument will be empirical and 
will consist in a rough quantitative correlation of magnetoresistance, 
Hall, and mobility data, based on (23). 

The curves of Fig. 8 show the first of the comparisons just mentioned. 
The abscissas are values of the ratio mw/mhq » where n H is the Hall mo- 
bility computed with (22) or (23) for a simple-model semiconductor, 
and ni, a is the corresponding quantity for pure acoustic scattering (6 = 
or 6 = 0). The upper two sets of curves show the values of the ratio of 
Hall to drift mobility and of the electronic Peltier heat relative to the 
band edge. The full curves were computed from (23), while the dashed 
ones were taken from the literature ' ' based on (22).* Note that 
the agreement is very good near mh/m«q = 1, but that, as one must expect, 
it becomes poor for mh/mhq < 2- For the magnetoconductivity (bottom 
curve) the agreement with the dashed curve, calculated from (22) by 
means of the tables of Beer, Armstrong and Greenberg, is even better. 
It is worth noting, incidentally, that Mansfield has also found mh/m to 
be insensitive to the form of the impurity-scattering law for n H /vHa > 
0.7 or so. 

Although it is encouraging to find such insensitivty of the predictions 
of impurity-scattering theory to the exact form of the scattering law, a 
comparison of the predictions of (23) with experiment is desirable as a 
test of its reliability. In applying (23) to n germanium we must, of course, 
take account of the anisotropy of the valleys, and of possible anisotropy 
in the impurity scattering. We shall do this in the same way as was de- 
scribed for the moment analysis of Section VI. Thus, we shall describe 

* Ref . 28 gives a general review and references to the earlier literature, Ref . 29 
an accurate table of integrals and Ref. 30 the most complete tables of Hall and 
magnetoresistance integrals. 
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Fig. 8 — Tests of the predictions of the impurity-scattering assumption (23). 
The full curves represent the predictions of (23), the dashed curves those of (22), 
for a semiconductor with isotropic effective mass. Abscissas are ratios of Hall 
mobility to its value for pure acoustic scattering. Ordiuates arc the ratios of the 
values computed for the designated quantities with impurity scattering to the 
values with pure acoustic scattering. The top curves pertain to n* — II, = (e 2 r)/ 
e(ei->, where n* is defined by (21); the acoustic value of n* - n e is 2kT/e. The 
middle curves pertain to m/z/V = <" 2 )(0/(eT) 2 ; the acoustic value is 3r/8. The 
bottom curves pertain to the low-field transverse magnetoconductivity ratio 
(c/m//#) 2 OWo-o) = (eT 3 )(eT)/(eT 2 ) 2 ; the acoustic value is 4/tt. The points represent 
empirical values for the various etched specimens as obtained from the observed 
ii///fi and the magnetoresistance constant b = (Ap/p# 2 )ioo 001 with abscissas based 
on the assumption; Mffa = 54.3, 38.6, 28.0, and 16.2 X 10 3 cm 2 /vs at 63°, 77°, 94°, 
and 131°K, respectively; both ordinates and abscissas represent ratios of observed 
quantities to values which would obtain if the u>\ , w* and w 3 of Table X were com- 
bined with the acoustic value of t± . The ratios of mh/m and of b to their ideal 
values are rather insensitive to the anisotropy of m*/r when this is allowed for in 
this manner. 
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the anisotropy by the w A , w* and W\ values given in Table X, and shall 
take the transverse magnetoconductivity of a [100] specimen to be given 
by 

Act 



^H/cY- (24) 

Ap 



[ (2 + ig a )(l +2w 3 )(2 + Wl ) l r (er J)(er*y \ 

L 3(1 + 2w 2 y- _ . (erJY J' 



pitxnH/cY 

an expression which follows from the known forms for A«r, <r and n H in 
the electron-group approximation, and the further assumption (19). 
Similarly, we shall take for the Hall and drift mobilities 

L (2 + «>i) _ _ (cTx) J 



/I a (2 + Wi) 



LTT 1 ' (2(5) 



so that the ratio of mh to its ideal acoustic-scattering value y. Ha can be 
obtained from the ratio of the right of (25) to its ideal value. If now we 
assume (23) for t ± , and the w n 's of Table X, we can use the curves of 
Fig. 8, since the ratio of the magnetoconductivity or the Hall mobility 
to the factor involving w„'s in (24) or (25) is just the value for the iso- 
tropic-mass case. 

As we have no direct measurement of m//o from which to compare (25) 
with (24) for a single specimen, we have adjusted this one parameter to 
lit the average behavior of the magnetoresistance of the purer specimens 
at 77.4°K. The mh« values for other temperatures were determined from 
the assumption m//.. a T -1 " 65 , this exponent having been chosen to make 
the fractional departure of n» from m«« scale a little more slowly than 
T~'\ as one expects from the Conwell-Weisskopf law. Using the ratio of 
the observed m// to this n, /a and the Wi and w 2 of Table X, we can deter- 
mine the ratio (eT±)(eT a± )/(eT ± )(eT a ± 2 ) of the last factor in (25) to the 
value it has when t x = r n ± , the value for pure acoustic scattering. With 
this ratio as abscissa, we have plotted, on the graphs of Fig. 8, the ratios 
of empirical value to acoustic value for the last factor in (24) and for 
the ratio of the last factors in (25) and (26), as obtained from the ratio 
of the low- and high-field Hall constants. Only data from etched speci- 
mens have been used, as extreme accuracy is required for a significant 
comparison; even so, there is a fair amount of scatter in the points. Much 
of this may be due to errors in m// resulting from errors — perhaps 1 per 
cent — in the assumed dimensions, or from the fluctuations in impurity 
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density; such errors will displace the \x a l\i points horizontally and the 
magnetoconductivity points along a line sloping at about 45° from up- 
per left to lower right. However, the curves account at least roughly for 
the at first sight surprisingly low magnetoconductivity of the pure speci- 
mens, and for the further reduction of magnetoconductivity and fx H /n 
found for the less pure specimen (Sample 601). 

In concluding this digression, a few words are in order regarding the 
validity of the assumption that nothing except ionized-impurity scat- 
tering interferes appreciably with the ideal acoustic scattering law r « 
t*. Optical-mode or intervalley scattering has sometimes been suggested 
as one of several possible causes of the departure of the mobility from 
a T~* law. Though such scattering freezes out at low temperatures, it 
is easily verified 13 that the amount of it which would be needed to account 
entirely for the mobility exponent could raise mh/m by 3 per cent or so 
at 131°K, depress it by perhaps 1 per cent at 77°, and lower the absolute 
value of the magnetoconductivity by 1 per cent or so. The fact that, for 
Sample 606 at 1P1°, the ratio n H /n is 0.992 of the acoustic-scattering 
value suggests that this and other departures from ideal acoustic-scat- 
tering behavior are small. 

7.3 Behavior of the Moments of H p \\,± 

Fig. 9 shows the way in which the various moments of n p , as defined 
by (20) for either component, would be affected by impurity scattering 



__i 



of the type (23), if n p were cc € ~' or e . Impurity scattering, of 
course, does not affect the zeroth moment, and if n p is independent of 
energy it does not affect any moment. We have chosen as abscissa the 
same ratio of Hall mobility to ideal acoustic mobility which was used in 
Fig. 8. Fig. 10 shows, analogously, the effect of impurity scattering on the 
second and third moments of the quantity (n* — n e ) « € ; the effect on 
the first moment has already been shown in Fig. 8. Comparison of Figs. 
9 and 10 with Fig. 8 shows that, for even our purest samples at 77°K, 
the effects of impurity scattering on the moments with n = 2 and 3 can 
be considerable if n p varies as rapidly as e~\ This is why we have un- 
dertaken such an elaborate discussion of these effects. 

The first step in analyzing the data of Fig. 7 is to correct n* — n e for 
the effects of impurity scattering and subtract it from the curves of Fig. 
7 to get curves representing the moments of H p || and U px , averaged 
over the purer specimens. To describe the average behavior of the purer 
specimens, we have chosen the values (cta )(er aJ .)/<€T X )(eT a j. ) = 0.89 
at 60°, 0.93 at 77.4°, 0.96 at 94° and 0.98 at 131°. These choices, though 
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Fig. 9 — Effect of impurity scattering on the moments n p (n) if n„ is assumed 
a t - " 2 (full curves) or « e - " 4 (dashed curves). Abscissas are the ratio of Hall 
mobility to its value for same anisotropy parameters Wi and w-> , but tj. <* e~ 112 . 

reasonable enough, differ slightly from other equally reasonable choices 
which one might make from Fig. 8; we shall discuss the sensitivity of the 
results to the choice in the next paragraph. From these abscissas in Fig. 

10 and the top curve of Fig. 8 we have determined the corrections to 
II* - n e (n> . The resulting values of n pJ . (n) and n p || (n> are shown in Fig. 

11 as full curves. The dashed curves are the moments to be expected for 
rip||. x « e°, e -I/4 , or e _1/2 , as determined from the values of im/nna just 
given and the curves of Fig. 9; the dashed curves have been fitted to the 
full curves at n = 1. 

A word about the sensitivity of the results to the choice of abscissas 
in Figs. 8 to 10 is in order. Shifting the abscissa 0.01 unit at 77° has a 
negligible effect on the moments with to = or 1. It shifts the full and 
dashed curves of Fig. 11 in the same direction for to > 1, the dashed 
curves being shifted much the greater amount; for to = 3, it shifts the 
e" 1 dashed curves by about 0.02 volt for n,,, , or 0.002 volt for U p± . 
At 131°K, the effect is more serious: a decrease of the assumed abscissa 
from 0.98 to 0.97, though affecting the to = 1 moments and n p || (3) very 
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Fig. 10 — Effect of impurity scattering on the moments n* — n,'" 1 . Abscissas 
are the ratio of Hall mobility to its ideal acoustic value for the same anisotropy 
parameters w\ and w> , but t± <x e _1/: . 

little, would lower the empirical — U p± {3) by about 0.001 volt, and would 
lower the corresponding ordinate of the e -3 curve for — n,,!!* 3 ' by 0.09 
volt. 

As is to be expected from our experience in constructing Fig. 7, the 
results for U P J are the more consistent. At all four temperatures the 
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empirical n pX curves lie between those for U p± <* e and c - ', and nearer 
to the former. By interpolation we estimate U p± « c - , the exponent 
being within 0.02 or so of this value at all four temperatures. The curves 
for n,,/"' have, for n > 1, roughly the behavior to be expected for 
n p| , cc r - 25 (77° and 94°K), e' or ° (60°) and e -050 (131°K); between 
n = and n = 1 the slope is steeper, corresponding on the average to, 
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Fig. 11 — Moments of n p j. and n pM (full curves), as determined from Fig. 7 
and the values of n* — 11, corrected for impurity scattering. The dashed curves 
show the behavior which would be expected if II^.j. were <x e", «-"« or e -1 ' 2 , with 
allowance for impurity scattering. 
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Table XI 
Computed fractional errors in the assumption that Q is the sum of the 
ideal Q e and Q p values which would occur in the absence of impurity 
scattering. Here B\Q\ = \Q |tme — I Q I ideal = 5 | Q e | + 5 | Q p | and 
5 I Qp I — &•» I Qp I + 5i | Q p | , where 5„, | Q„ | arises from the departure 
of Wy = m* x (eT l \)/m*i\(eT ± ) from its acoustic-scattering value, and 
5i | Q p | arises from the corresponding departures of the lip/", IT, 
denned by (20). 



(i) 



Sample 


Estimated 
Abscissa for 
Figs. 8 & 9 


T, °K 


*I0.|. 

nv/deg. 


Sw I Q P 1 , 

j»v/deg. 


Si 1 Q P i , 
Mv/deg. 


5IQI 

IQI 


606 


0.91 


60° 


+5 


^+5 


-25 


~ -0.005 




0.99 


122° 











0.000 


603 


0.93 


77° 


+3 


^+5 


-4 


< +0.002 


601 


0.75 


94° 


+ 18 


+ 16 


-17 


+0.010 


596 


0.56 


87° 


+38 


+57 


-30 


+0.041 



say, e ° 4 at the three lower temperatures. Considering these variations 
and the random scatter of the points in Fig. 7, we can state as our best 



estimate that n p || « c " 



at low energies and that its energy de- 



pendence probably becomes more rapid at higher energies. There is no 
reliable evidence that the energy dependence of either component 
changes appreciably in the temperature range 60° to 131 °K, although 
the accuracy of the determination of energy dependence is not very good 
for IIp|| . The ratio nj,|| (1) /n p x (1> seems to be decreasing slightly with in- 
creasing temperature, although here again it is hard to be sure that the 
effect is real. 

7.4 Effect of Impurity Scattering on the A ssumed Q 

As we have mentioned in Section V, all our thermomagnetic quantities 
were measured relative to the value of Q at H = 0, and the values orig- 
inally tabulated for them were based on the value of Q characteristic of 
the observed carrier concentration and pure acoustic scattering. If im- 
purity scattering causes the total Q to depart appreciably from this as- 
sumed value, these values will require correction. The results obtained 
in this section allow us to estimate how much effect impurity scattering 
has on Q e and Q p , and so to check this point. 

According to Table IX, impurity scattering can affect Q p in two ways: 
through changing the anisotropy w(=Wi) of m*-T~ , and through its 
effect on the moments n p || (1) , n pJ . (1) . The effect on Q e is, of course, simply 
obtained from the top curve of Fig. 8. Table XI shows some sample eval- 
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uations of the two kinds of effects on Q p and of the effect on Q e , computed 
from Figs. 8 and 9, the w's of Table X, and the assumptions n p || a 
e -0 ' 30 , U p± a e -08 . It will be noted that the change in w largely compen- 
sates the changes in n p ||, x (1> ; this explains the previously puzzling lack 
of sensitivity of Q p to impurity scattering. The fractional errors 6Q/Q in 
our assumed Q values are seen to be at most a fraction of a per cent for 
the purer specimens ; this implies errors in the ordinate of Fig. 7 of rather 
less than a per cent. For Samples 601 and 596, however, at the tempera- 
tures shown, 8Q/Q is 1 per cent and 4 per cent, respectively, and 511/ 
(n* — II) is about twice as great. 

7.5 Confirmation of Impurity-Scattering Theory on Sample 601 

We have based our corrections for impurity scattering in the purer 
samples on plausible theoretical reasoning and on favorable though not 
ideally consistent evidence from isothermal electrical measurements. It 
is therefore of some interest to see if the thermomagnetic measurements 
on a sample with a sizable amount of impurity scattering differ from 
those on the pure samples, in the way in which we would predict on the 
basis of (23) and Figs. 8 to 10, assuming the energy dependences of 
rip|| and n pX which we inferred from Fig. 1 1 . 

Fig. 12 shows the comparison just mentioned for Sample 601 at 94°K. 
This sample was chosen for the comparison because good isothermal data 
were available for it and because it had about the maximum amount of 
impurity scattering for which the predictions of (23) are likely to be re- 
liable. (For some of the other samples, the electrical and thermomagnetic 
data were at different temperatures, and for 596 at 87° the divergence 
between the full and dashed curves of Fig. 7 is beginning to be appreci- 
able.) The dashed curves in Fig. 12 are the curves of Fig. 7, and repre- 
sent the mean moments of the total n* — YL\.± for the purer samples. 
The plotted points represent the empirical moments for Sample 601, com- 
puted from the data of Table VI in the same way as for the pure sam- 
ples of Fig. 7. The small temperature-gradient correction just discussed 
is, of course, included in the corrections of Table VI. These points are 
to be compared with the full curves, which were constructed by comput- 
ing the differences n l± (n) (601) - Yl\\,± ln) (ideal) theoretically and 
combining them with the dashed curves. The differences in question 
were computed by the following steps: (a) From (25), the observed 
Hh of 601, the empirical W\ and w 2 of Table X, and the assumed acoustic 
value n„ a (94°) = 28,200 cm 2 /vs (probably good to within a few per 
cent), we computed the abscissa appropriate to Sample 601 in Figs. 8 to 
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Fig. 12 — Comparison of the empirical moments II* — Hm.j.'' for Sample 601 
(points) with those of the purer samples of Fig. 7 (dashed curves) at 94 °K. The full 
curves are moments computed from the latter, with allowance for the amount of 
impurity scattering present for Sample 601. 



10. (b) From these figures and the assumptions II ^ <x c 030 t U p± ex 
e" 008 , the ratios of n ph J n) (601) to n ph J n) (ideal) were obtained, (c) The 
same ratios were obtained for the average of the pure specimens, as- 
suming an abscissa of 0.97 in the figures. (This is larger than corresponds 
to the above /*//« ond the mobilities measured on the unetched samples, 
but is more reasonable in the light of the 77° data and the effect of etch- 
ing at that temperature.) (d) From steps (b) and (c) and the empirical 
n P ||. J . ( " ) of Fig. 11 the n P ||,j. (n) to be expected for Sample 601 were cal- 
culated, (e) The theoretical U* — Tl e n were calculated for Sample 601 
from Figs. 8 and 10, and combined with the n p | lX just obtained. 

It will be seen that, for Hjj C especially, the empirical points fall 
much closer to the curve computed with allowance for impurity scatter- 
ing than to the curve for the pure samples. Since the theoretical differ- 
ence between the two curves is proportional to the rate of variation of 
IX|| with energy and to the departures of the ordinates of Figs. 8 to 10 
from unity, the agreement is a significant bulwark for the analysis of this 
section. However, it is of more qualitative than quantitative significance, 
for the inputs to the calculation were numerous and some of them a little 
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uncertain. The discrepancy for TT /(J . (3) may well be due to the great sen- 
sitivity of the II* — II,. <3> curve of Fig. 10 to the abscissa chosen. 

VIII. FURTHER OBSERVATIONS SUBSTANTIATING AND EXTENDING THE 
PRESENT MODEL 

In this section we shall discuss several parts of our data which are not 
adapted to as precise an analysis as that of Sections VI and VII, but 
which can still be semiquantitatively explained and give indications of 
how well our conclusions can be extended to higher energies and higher 
temperatures. 

8.1 The High-Field Asymptote of BH 

So far, we have made no use of our observations of the Nernst coefficient 
at large magnetic fields. If the asymptotic high-field behavior were ac- 
curately known, it could provide a valuable extension of the moment 
analysis of Section VI. For, as (80) of Appendix B shows, this asymptotic 
behavior involves the moments (20) for n = and — 1 , so that measure- 
ments in two orientations could be used to determine n/" 1 ' and n.L (-1) . 
In practice, however, the accuracy of such a determination would be 
very questionable, as it would depend on the already rather uncertain 
zeroth moments and on accurate knowledge of the asymptotic BH, 
which, as we shall presently see, is rather uncertain for various experi- 
mental reasons. We have therefore chosen merely to compare the high- 
field behavior of BH with the asymptotic behavior predicted by the 
special -case formulas of Table II, with the object of showing that there 
is at least fair agreement with the model previously deduced. 

As we have noted in Section V, the small value of the Nernst constant at 
high fields makes it very sensitive to falsification by the Hall fields of the 
longitudinal counter-currents which will be set up if surface conduction 
partially short-circuits the thermoelectric field; a surface conductance of 
a fraction of one per cent can be serious. We did not realize this soon 
enough to carry out extensive Nernst measurements on etched specimens. 
The lowest temperature at which such measurements were made was 
94°K. At this temperature the difference between the behavior of Sample 
(506 etched and unetched was profound, the unetched having about twice 
the Nernst field of the etched at 18 kilogauss. At higher temperatures, 
the surface effect gets even worse, because the counter-currents are driven 
by the full Q, including the term proportional to the distance of the Fermi 
level below the band edge, whereas the Nernst effect arises only from 
Q p and the small transport term in Q e . 
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Fig. 13 — Comparison of the asymptotic behavior of the Nernst field constant, 
BH, at high fields with the predicted (B e -f- B P )H of the simplified asymptotic 
formulas of Table II. The three experimental curves are for two samples of differ- 
ent orientations; note that the two scales are different. The short-dashed curves 
are computed from the formulas of Table III for U p independent of energy, the 
long-dashed curves for H p <* « _1/2 . For both sets, the values w = 0.061, p = 9.8 were 
assumed. The unit of BH is volts/degree. 



Fig. 13 shows a comparison of three observed curves of BH against H 
with the leading (If 1 ) term in the asymptotic expression calculable from 
Table II. The data for Sample 606 at 94°K are the most reliable, since 
this sample has been etched; however, these data do not extend to as 
high a value of nH/c as is obtainable at lower temperatures, so there 
may be some doubt as to how well they approach asymptotic behavior 
at 18 kilogauss. The data for Sample 606 at 60°, though on an unetched 
specimen, appear not to have been very seriously falsified by surface 
conduction, since they match the 94° data pretty well — as far as they 
go — if ordinates are scaled proportionally to Q p and abscissas to ju, 
and if a little allowance is made for the electronic contribution. However, 
at the high-field end, the 60° data lie closer to the theoretical asymptote 
for n p cc c~' than to that for n p independent of energy, whereas the 
reverse is the case for the 94° data; this may be due to a slight raising 
of the experimental curve by surface conduction. 

The data for Sample 603 — with H in a [1 10] direction — show the 



V 
PHONON-DRAG EFFECT IN n-TYPE GERMANIUM 715 

vastly slower saturation which is to be expected from the fact that, for 
this orientation, half the valleys have a very large cyclotron mass. Since 
the curves for the two samples are plotted on different scales, the values 
of BH at 18 kilogauss and 60°-61°K actually differ by a factor of 14 in 
the two orientations, roughly the factor predicted by the theory. As 
nearly as one can estimate from the range available, the experimental 
curve seems to be heading for an asymptote closer to that for II P inde- 
pendent of energy than to that for II P cc e ~ 2 — perhaps even the other 
side of it. This is a slight discrepancy with theory which might again be 
due to the presence of a little surface conduction. 

Thus, it appears that the qualitative features of the high-field Nernst 
curves are accounted for by the theory. The curves also give some in- 
formation on the value of the anisotropy p = n P ||/n pJ . . Since the en- 
tries in the next to the last column of Table II are almost proportional 
to (p — 1), it is probably safe to say that the present curves indicate 
a p within 20 per cent or so of the value derived in the last section. Re- 
garding the variation of the weighted average of II P || and U p± with en- 
ergy, about all that can be said is that this variation again appears not 
to be as rapid a decrease as e - . 

8.2 Reversal of the Sign of B with H 

All our analysis so far has been based on data at temperatures in the 
range 60° to 131°K. No comparably accurate analysis can be made 
at higher temperatures, partly because the limits of the various quanti- 
ties as H — > oo can no longer be determined, and partly because the 
value of Q p is much less accurately known than in the lower range. What 
evidence there is, however, suggests that, at temperatures in the range 
1 50° to 235°K, the anisotropy and energy dependence of n p are similar 
to those which we have found in the liquid-air range. This evidence is 
provided partly by the magnitudes of the low-field B and AQ/H , which 
we shall discuss presently, and partly by an interesting reversal of the 
sign of B with increasing magnetic field. The latter effect occurs at tem- 
peratures at which the negative B e (0) and the positive B p (0) almost 
cancel each other. This cancellation occurs near 175°K, B(0) being 
greater than zero below this temperature, less than zero above it (al- 
though easily falsified by surface conduction). Fig. 14 shows some curves 
of B versus H at temperatures in this range, taken on etched samples. 
The data previously reported 1 on 606 and 604 unetched showed no region 
of negative B(0), but now both show sign reversals at 171° to 174°. 
What we wish to discuss now, however, is the fact that, as H increases, 
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Fig. 14 — Variation of the Nernst coefficient B with magnetic field H, for two 
orientations, showing the sign reversal which takes place at temperatures near 
where B e and B p compensate each other. For comparison, the computed 5.(178°) 
is —4.0 X 10 _9 v/deg gauss. 

the B of Sample 604 goes from negative to positive, and that of 606 
from positive to negative. 

Although we have not elaborated the theory of B at intermediate 
fields, it is doubtless safe to assume that, for all field strengths, B p > 
and B e < 0. We can get a fair idea how the ratio of B p to B e varies 
with H by comparing the low- and high-field limits of | B p /B e | as given 
by the formulas of Tables I and II. Table XII gives the comparison, for 
several sets of assumptions about the anisotropy and energy dependence 
of lip . We have constructed schematic curves of B e and B p against H, 
shown in Fig. 15, from this information and the knowledge 1 that B p 
decreases monotonically with H for H || [001], while it has an initial 
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Table XII 

Comparison of high- and low-field ratios of B p (positive) to B e (nega- 
tive), for several orientations and for different assumptions regarding 
the anisotropy p = n p ||/n p j. and the dependence of n p on energy e. 
The numerical values given are for iv = m* X Tf/m*Tu ± = 0.06. 



Direction of 



[001] 
1011] 

[Oil] 



\T 



J_H 

[100] 

[oil] 



Ratioo Mi^)r o lral ,for 



n„ independent of «, any p 



i** - 6.2 

3w 

(2 + w)(l + 2m*) 
3mi(1 + 7w + w 2 ) 



= 8.1 



P = s 



0.45 
1.46 

2.27 



P = 10 



0.50 
2.06 

3.31 



p = 20 



0.54 
2.58 

4.17 



rise for H || [011], if p = n p ||/n p j. is large. From these curves it is clear 
that, if the anisotropy and energy dependence of n p are similar at the 
temperatures of Fig. 14 to what we found in Fig. 11, then B should be- 
come more negative with increasing H for Sample 606, more positive for 
604. This is exactly what is observed. 

Only the crudest sort of quantitative limits can be derived, from this 
comparison, for p and the rate of variation of n p with energy. From an 
inspection of Table XII, Fig. 15 and the curves of Fig. 9 of Ref. 1, it 
seems clear that the sizes of the shifts of B with field in Fig. 14 are 
eminently consistent with the previously found anisotropy and energy 
dependence, and that they would be hard to reconcile with a n p that 
increased with increasing energy, or with a p close to unity. Specifically, 
Table XII suggests that the theoretical curve for 606 is nearly inde- 
pendent of p, so that the downward slopes in Fig. 14 must arise from an 
inverse energy dependence of n p . The swing in B between zero and 
5000 gauss is of the order of a twentieth of | B e (0) |; this would not be 
unreasonable for n p « e~*. The swing for 603, on the other hand, is 
probably due to the fact that B p rises with H while B e falls; if so, it 
should be of the order of one or two times the rise in B p , a rise due 
mainly to the anisotropy of n p (p » 1). For the observed anisotropy at 
liquid-air temperature (p « 10) this rise (see Fig. 9 of Ref. 1), at the 
HhH/c corresponding to 5000 gauss at 178°, is perhaps 15 per cent of 
B p (0). The observed swing is about twice this. 
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Fig. 15 — Schematic variation of B. and B p with magnetic field, for two orient- 
ations, assuming the temperature such that B e (0) = — B p (0). Two B p curves are 
shown for each case, corresponding to n p <* e » and e _m , respectively. It is assumed 

that n P n/n p j. » l. 

8.3 Analysis of High-Temperature Data 

The sort of moment analysis which we made in Figs. 7 and 11 gets 
increasingly difficult as T increases, for several reasons: high-field limits 
become unavailable; the phonon-drag effects become a much smaller 
part of the total; because of this, it is easier for such things as surface 
inversion layers to falsify measurements of bulk thermomagnetic effects. 
We shall nevertheless try to draw what inferences we can from the tem- 
perature at which B(0) changes sign, and from the totality of measure- 
ments on Sample 601 at 234°K, a sample pure enough to be fairly free of 
impurity scattering at this temperature, yet sufficiently highly conducting 
to be insensitive to surface conduction. 

As Fig. 14 shows, B(0) reverses sign in the pure, etched specimens at 
174°K (603) or 172°K (606). From the average of these and interpola- 
tions of the Q p values of Table VI we can deduce a value for the f p of 
(13). With the acoustic-scattering value of B e , we obtain f p = 0.226, 
a value close to the average f p of 0.229 found for the pure specimens at 
77°K, corrected for the effect of impurity scattering on B e . This suggests 
either that the anisotropy and energy variation of TL P are very similar at 
173° and 77° or that the difference in anisotropy just happens to be 
such as to compensate for the difference in energy dependence. The 131° 
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Fig. 16 — Results of an attempt, to determine empirically the moments Up\\,x m 
defined by (20) for Sample 601 at 234°K. The points O, + and X correspond re- 
spectively to three choices of the assumed n p ± (1) ; it will be seen that the O choice 
gives the most nearly plausible behavior for n P n (n) . 

data of Fig. 11 suggested a possible decrease of the anisotropy with 
increasing temperature, and a possible increase in the rapidity of varia- 
tion of IIp|| with energy. These factors would affect f p in opposite ways, 
so they could cancel. 

We have attempted a moment analysis like that of Sections VI and 
VII for Sample 601 at 234°K. We adopted the plausible anisotropy 
Wl = Wz = Ws = 0.056, (see Table X and the discussion of Section VI, 
especially the remarks about change of w with temperature). This is a 
value which happens to fit the anisotropy of the magnetoresistance. The 
data give II* - U ± (3) = 0.0507 volt, n* - II|, (3) = 0.298 volt. The only 
additional knowledge obtainable, however, is one relation between H|| ' 



and n x (1) (to fit Q) and one relation between n,, and n x l (to fit B] 



(n) 



If, as before, we assume II /"' to be linear in n, we can adjust II x until 
the curve of II|| (n) against n is reasonably smooth. Fig. 16 shows the 
results for n„||, J . (n) after subtracting off the values of n* - n e (n) , slightly 
corrected for impurity scattering. Several choices of n pJ . (1) have been 
made, and their consequences are shown in the figure. There are strong 
suggestions that n^^/IIp., ' is still smaller than at 131 °K, that n p| | 
varies more rapidly with energy than at 77° and that U p± continues to 
vary slowly. However, the data are so limited that only the last of these 
points can be considered reasonably established. 



IX. SUMMARY AND IMPLICATIONS OF THE RESULTS 

9.1 Evidence on n p y(e) and II pX (e) 

We have concluded that, near liquid-air temperature, the dependence 
of the phonon-drag Peltier tensor of an energy shell on the energy e of 
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the shell is given, at least roughly, by n p ||(e) « e -0 ' 3 , n pJ .(e) <x e~ ' , 
while the first-moment ratio <€ , n p ||)/(e Yl p ±) is about 9.6. These con- 
clusions were based primarily on the analysis of high- and low-field AQ 
data and low-field Nernst data on the several purest samples, with aux- 
iliary use of mobility and magnetoresistance data. The use made of the 
low-field data depended to some extent on corrections for impurity 
scattering, the method of making the corrections having been roughly 
checked both against purely electrical measurements and against alter- 
native theoretical assumptions. The data for Sample 601, which showed 
a sizable effect of impurity scattering, were consistent with the assumed 
characteristics of n p (e) combined with the assumed impurity-scattering 
law. The conclusions regarding n pJ .(e) were very insensitive to the more 
uncertain of the inputs and were consistent from sample to sample and 
from one temperature to another. Those regarding n p ||(e) were more 
sensitive and showed more fluctuation, but always showed roughly the 
same behavior. High-field Nernst data, though less suited to accurate 
analysis, independently indicated the ratio of II P || to U px to be of the 
order of 10, and also indicated a decrease of the weighted average of these 
quantities with e, slower than e . 

The high-field limit of the transverse AQ's gave strong but not entirely 
certain indications that the decrease of n P || with increasing energy be- 
comes more rapid than the above-mentioned e~ 03 at high energies. The 
variation of the anisotropy and energy dependence of n p with tempera- 
ture was too small to be detected over the range 60° to 94°K. How- 
ever, we found some rather uncertain indications that n p ||/n pJ . decreases 
slowly with increasing T, and that the energy dependence of n p || be- 
comes steeper. These conclusions came from a comparison of 131° 
thermomagnetic data with data at lower temperatures, and from studies 
of low-field thermomagnetic data at 234°. The low- and intermediate-field 
behavior of the Nernst coefficient near 175°K was found to be consist- 
ent with the anisotropy and energy dependences found from the previous 
sources, but did not suffice to check the temperature variation of these 
two factors individually. The data at all temperatures are clearly in- 
consistent with a law n p || <* e 5 , such as would result (see below) if 
the relaxation times of the phonons were independent of frequency. 

9.2 Implications for Phonon-Phonon Scattering 

The energy dependence of n p (e) reflects the dependence of the relaxa- 
tion times T p h of the different phonon modes on their frequency o>. For 
example, if all branches of the phonon spectrum had 

Tph = w' X (function of direction), (27) 



PHONON-DRAG EFFECT IN n-TYPE GERMANIUM 721 

then, since the wave numbers, or frequencies, of the modes coupled to 
electrons of energy e scale proportionally to e\ we would have, as « 
varies 

n„ <* r _1 X (a mean t p u) 

J+*« (28) 



a e 



If the longitudinal and transverse branches were to obey (27) indi- 
vidually, but with different exponents h and t t , then the effective energy 
exponents for Yl Pll and H px would be intermediate between ^ + \tt and 
| + %i t , and the exponents would in general be different for n p || and 

n pi . 

Our liquid-air-temperature result is thus that the phonons interacting 
with a current parallel to the valley axis have an average relaxation time 
T P h « w -1 '", while those interacting with a current in the perpendicular 
direction have an average r p & « aT 1 ' 15 . These are to be compared with 
the "ideal" phonon-phonon scattering laws ' " 

longitudinal branch: t p i, « co ", (29) 

transverse branch: t p u <* to , (30) 

valid if co « the frequency of the average thermal phonon and if only 
threc-phonon noncollinear collisions contribute to the relaxation of the 
modes concerned. Estimates 32 of the frequency at which appreciable de- 
partures from (29) or (30) set in have suggested that, for germanium at 
liquid-air temperature, these laws should be approximately obeyed, but 
that departures from them may well be appreciable. As it can be shown 
that T ph goes as a higher power of co when co gets large, it seems likely 
that the finiteness of co will lead to higher, rather than lower, negative 
exponents. Thus, at present a more plausible explanation of the observed 
exponents is that the longitudinal and transverse branches contribute 
comparably to the phonon-drag phenomenon. Since it is known' that the 
longitudinal and transverse branches contribute comparable amounts to 
the scattering of the electrons, this would require that, in the frequency 
range most important for the scattering of electrons, the relaxation times 
of the transverse modes be comparable with or rather larger than those 
of the longitudinal modes. This is contrary to what one would anticipate 
from (29) and (30), but cannot yet be considered unreasonable, in view 
of our present ignorance about scattering matrix elements, etc. 

There are two further possibilities which would lead to departures 
from the asymptotic laws (29) and (30). One of these is three-phonon 
collisions in which the three phonons belong to the same branch and have 
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almost collinear wave vectors. If such collisions are possible, they will 
change (29) into the same form as (30). For this to occur for a low-fre- 
quency longitudinal mode, it is necessary that some higher-frequency 
mode with the same direction of propagation have a higher group 
velocity than a mode of infinitesimal frequency. 32 This could happen, for 
example, if the curve of frequency against wave number in the longi- 
tudinal branch curved upward instead of downward. Although this does 
not seem impossible a priori, no real crystal or mathematical model has 
to our knowledge been found with this property. Moreover, by now our 
knowledge of the elastic spectrum of germanium 33 is good enough to 
make it fairly certain that this does not occur for the present material. 

More intriguing is the possibility of higher-order processes, i.e., colli- 
sions involving four or more phonons, 34 or possibly even processes more 
appropriately describable as relaxation effects. 35 From the standpoint of 
collision theory, a calculation of, say, three-phonon collisions becomes 
invalid when the frequency of one of the modes involved becomes less 
than the natural widths t p a _1 of the other modes, or when the failure of 
energy conservation in an "energetically impossible" process gets less 
than some relevant hr^T 1 . Moreover, even when neither of these catas- 
trophies occurs, and the three-phonon process that is being considered is 
perfectly meaningful, it may happen that the number of such possible 
three-phonon processes is so small that the relaxation of the low-fre- 
quency mode involved occurs predominantly by higher-order processes. 

However, the rates of four-phonon processes vary much more rapidly 
with temperature than do those of three-phonon processes. At the fre- 
quencies of the phonon-drag modes, which are high compared to the re- 
ciprocal thermal relaxation time of the majority modes, relaxation can 
be described as due to some, though not all, of the collisions involving 
four or more protons. Since the temperature dependence of Q p in the 
range 60° to 90°K is not far from that to be expected for three-phonon 
processes, 2 we must conclude that higher-order processes are at most a 
small correction in this range. But if they are at all appreciable at one 
temperature, they must become completely dominant at a higher tem- 
perature. The changeover from dominance of three-phonon processes to 
dominance of higher-order processes would almost certainly change the 
anisotropy or the energy dependence of n p ( e) . In particular, if the relaxa- 
tion time of the phonon-drag modes became independent of frequency — 
as it would for the simplest type of higher-order processes — n p would go 
as «*, in violent disagreement with our observations at all temperatures. 
The absence of any abrupt change of the energy dependence or aniso- 
tropy of n p with temperature thus shows that, for the phonon-drag 
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modes, three-phonon processes are more important than higher-order 
processes, up to at least 234°K, and that over most of this range the 
dominance of the former must be extreme. 

It is not possible to follow the phonon-phonon scattering to tempera- 
tures much below 60°K without a theory that takes account of boundary 
scattering, an effect which spoils the crystalline symmetry of the prob- 
lem. In fact, some of the apparent changes which Fig. 11 shows in the 
behavior of n P || at 60° may well be due to the fact that boundary scatter- 
ing is beginning to be too large to be corrected for by the crude method 
of Section V. 

9.3 Transport Theory and Deformation-Potential Theory 

Although the study of the purely electrical properties of n-type ger- 
manium was only an incidental feature of the present work, our results 
have led to a number of new conclusions in this field, which we shall 
summarize here. 

To begin with, we have found considerable evidence (presented in 
Appendix A and summarized in Section II) that in a magnetic field one 
need not abandon conventional transport theory in crystal-momentum 
space in favor of a quantum treatment until the cyclotron spacing Tito c 
has become close to kT. Specifically, we found that, at Tiu) c /kT ~ ^ and 
for acoustic scattering, the longitudinal and transverse resistivities are 
correctly given by conventional theory to within 5 per cent or so, and 
in most cases probably rather closer than this. Even at fico c /kT tt 1 our 
AQ evidence suggests, again for predominantly acoustic scattering, that 
the errors may be only of the order of 10 per cent or so. 

We have found the magnetoresistance and Hall data at high and low 
magnetic fields to be fairly consistent with the predictions of the electron- 
group theory, i.e., with the picture which approximates the effect of the 
scattering processes by relaxation times r\\(e) and Tj.(e). For the most 
accurately measured cases (etched samples at 77.4°K), the accuracy of 
the fit can be described as follows: The single parameter Wi = m ± *(eT\\)/ 
m\*(eT ± ) can be chosen to reproduce the longitudinal high-field mag- 
netoresistance limits p( co )/p(0) for the four purest samples in [100], 
[110] and [111] orientations to within the reproducibility of measure- 
ments from sample to sample (2 per cent or so). The values of p( °o )/p(0) 
for three transverse orientations also agree with predictions based on 
this anisotropy and the assumption t\\,± <* « , to within a somewhat 
larger uncertainty (worst discrepancy 13 per cent), due in part to sur- 
face conduction, etc., and in part to quantization. The low-field magneto- 
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resistance constants obey the symmetry relation characteristic of the 
electron-group approximation (the relations at the bottom of Table IV 
with Ap substituted for AQ P ) about as accurately as they obey the re- 
quirements of macroscopic symmetry (consistency to a few per cent). 

Concerning the energy dependence of ry and t± , the most accurate 
evidence we have is that from the ratio of low- to high-field Hall constant, 
i.e., the ratio of Hall to drift mobility. As was shown in Section VII and 
Fig. 8, this ratio, when measured on etched samples, gets closer and 
closer to the value predicted for r^j. « «~* as the importance of im- 
purity scattering is reduced. We were able to conclude that only a very 
small amount of optical-mode or intervalley scattering, or of any other 
departure from the law t|,j. « e~*, can be present at temperatures in 
the range 77° to 131 °K. 

In the calculation of the absolute value of the acoustic mobility, n„ , of 
n-type germanium from deformation-potential theory, 3 two deformation- 
potential constants are required. One of these can be obtained with some 
accuracy from piezoresistance measurements, and the other can be in- 
ferred if magnetoresistance or other data give a reliable value of the 
acoustic-scattering anisotropy tj/t x . Our present measurements have 
yielded much better values, both for n a and for w = m ± *T\\/m\\*Tj. , 
than were available when Ref. 3 was written, so it is of interest to see 
how well the deformation-potential calculation of y. a now works out. 

The first step is to obtain t\\/t ± from the cyclotron masses 36 m*.* = 
0.082m, m y * = 1.58m (assumed for the moment to be valid near liquid- 
air temperature) and the high-purity limit w = 0.061 , presumably 
characteristic of acoustic scattering. The result is t||/t x = 1.18. From 
this and Equations (52) and (53) or Fig. 5 of Ref. 3 (also based on the 
masses mentioned), we get the value —0.40 for the ratio Ed/E* of the 
deformation-potential constant for stretch along a direction normal to 
the valley axis to that for the combination of a compression along such a 
direction with an equal stretch along the valley axis. From this value of 
Sd/H,, and the piezoresistance constant " 

-'"' = !!m = 30 ' 400/r (31) 

we obtain Ed , S« , and thence the acoustic mobility. The value pre- 
dicted in this way for 77.4°K is 

n a (77.4°) = 45,100 cm 2 /vs. (32) 

The empirical value is obtained by dividing the ideal value liHa/v-a = 
0.933 into the observed Hall mobility, after correction for impurity 
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scattering by the estimates described in Section VII. With the value 
38,600 cm /vs there adopted for n Ha (see Fig. 8), we get 

n a (77.4°) = 41,400 cm 2 /vs. (33) 

The agreement between (32) and (33) is gratifying. Although the 8 
or 9 per cent discrepancy between them is larger than the uncertainty 
in either due to uncertainties in the experimental inputs (perhaps 2 or 
3 per cent), it is of the right sort to be attributable to departures of the 
masses or deformation-potential constants at 77° from the low-tempera- 
ture values we have assumed; moreover, the inherent inaccuracy of the 
deformation-potential method could conceivably be of this same order. 
As regards the former point, we may notice that the difference between 
mobilities varying as T~ ' and as T~ u 5 amounts to 8 per cent between 
46° and 77°K. If the latter variation, which we found to be closer to the 
truth over the range 00° to 94°K, were due to something like a tempera- 
ture variation of the masses which disappeared below 40° or so, the dis- 
crepancy would be eliminated. 

Regarding impurity scattering, finally, we have found in Section VII 
that it is possible to correlate with each other the effects of ionized- 
impurity scattering on mobility, Hall-to-drift-mobility ratio, magneto- 
resistance, etc. — at least in the range where impurity scattering is small 
to moderate. In this correlation the anisotropy of the impurity scattering 
was treated as an empirical parameter; beyond this, no specific assump- 
tion was necessary about the mechanism of impurity scattering, except 
that it affected low-energy carriers much more than high-energy ones. 
The correlations we were able to make in Fig. 8 and Fig. 12 seem to be 
quantitatively significant, but still leave something to be desired. 

All these items deserve more careful study. Our Table XV (Appendix 
B), incidentally, which gives Maxwellian averages of powers of the 
energy, may prove handy for various types of calculations. 

9.4 Surface Effects 

We have encountered serious effects of surface conduction in a number 
of places in the present study. Our experience justifies the statement that, 
at all temperatures, one must eliminate or allow for surface effects before 
one can get really accurate measurements of any bulk electrical property 
of high-purity material. For some properties, even the qualitative be- 
havior can be completely falsified by surface effects. Specifically, we have 
found that a surface contribution of less than one per cent to the con- 
ductance can greatly affect the apparent high-field transverse magneto- 
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resistance and the high-field Nernst coefficient near liquid-air tempera- 
ture and below, and that it can change the sign of the low-field Nernst 
coefficient in the range 175° to 250°K. Near liquid-air temperature, 
changes in surface treatment can alter the apparent conductivity and 
Hall coefficient by amounts ranging up to several per cent. 
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APPENDIX A. ORBITAL QUANTIZATION 

We give here the reasoning behind a number of statements which one 
can make, including those made in the text. We shall take up first some 
quantitative statements about the relation of the magnetoresistance of a 
multivalley semiconductor to that of the simple model assumed in the 
existing theoretical literature (items 2 and 3 below). Next we shall pre- 
sent empirical and theoretical evidence on the smallness of the effect of 
quantization on longitudinal magnetoresistance (item 4). Finally, in 
items 5, 6 and 7, we shall discuss the relation of the effect of quantization 
on AQ to that on magnetoresistance, and shall present evidence that 
these effects are fairly small even for transverse orientations. As the 
basis of the whole treatment we shall take the fairly obvious statement: 

1. For the region iiH/c » 1, the orbital-quantization correction functions 
for p(H)/p(0) and Q p (H)/Q p (0) are practically the same functions of H 
as they would be if the electron-phonon coupling were allowed to approach 
zero. 

In other words, the difference between the true p(H)/p(0) or Q P (H)/ 
Q p (0) and the value calculated without quantization is a function of the 
cyclotron spacing hco c (primarily of fuo c /kT), which depends, of course, 
on the anisotropy of the electron-phonon interaction, but which depends 
on the magnitude of the latter (i.e., on ixH/c) in such a way as to ap- 
proach a limit uniformly as the scattering goes to zero (nH/c — > <»). 
This is, of course, what one expects intuitively; it follows from any trans- 
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port theory in which p(0), p(H), etc. are representable by expansions in 
powers of the strength of the scattering. 

We shall show first that, if the scattering processes are energy-conserv- 
ing collisions with a squared matrix element that is the same for all 
transitions within a valley, the following statement holds: 

2. In the approximation of isotropic scattering over an energy shell 
(i.e., T|| = Tj.), the longitudinal magnetoresistance ratio p(H)/p(0) for 
any multivalley case equals the product of the value calculated in the ab- 
sence of quantization by a weighted average of the values calculated for iso- 
tropic m*'s equal to the cyclotron masses of the various valleys for the given 
direction of H, provided y-H/c ^> a value which for germanium is of the 
order of a fraction of m\*/mj.*. 

Thus, for longitudinal magnetoresistance, the orbital-quantization cor- 
rection can be deduced from the simple-model calculations in the litera- 

. ^ 6,6,7,8 

ture. 

For the proof we note first that when E || H in, say, the z direction, 
p zz = l/<r Z2 ; this is exact for all H in symmetry directions, and asymp- 
totically exact as // — > co in any direction. Therefore we need consider 
only <r« , which is a sum of contributions <r« (,) from the different valleys 
i. For any such valley we can reduce the transport problem to that for 
the simple model by means of a transformation which projects the energy 
surfaces into spheres. 3 In considering this transformation we must re- 
member that it is not an orthogonal one, and that it therefore carries 
with it a need to distinguish between covariant and contravariant vec- 
tors. For covariant vectors, F, the transformed vector, F, is related to 
the original vector F by 

F|| = r% , F A = F x , (34) 

where || and _|_ refer to the valley axis and r = m ± */m{\*. For contra- 
variant vectors G the relation is 

On = r*G,| , G, = G. (35) 

The scalar product of a covariant and a contravariant vector is un- 
changed by the transformation. It is easily verified that, if we regard E, 
the crystal momentum P, and the vector-potential A as covariant vec- 
tors, while regarding velocity (or j) and H as contravariant vectors, all 
the usual relations involved in the Schrodinger equation and the trans- 
port equation are the same for the transformed and original vectors, 
except that in the transformed space the effective-mass tensor becomes 
isotropic. 

Now let aiCH.) be the conductivity in the transformed system when 
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E || H, and a- t (K) be that for E _L H. These conductivities will be inde- 
pendent of the orientation of H or of the azimuth of E about H, if the 
scattering for // = is isotropic over each energy shell, as is known to 
be very nearly the case for acoustic scattering in n-type germanium, and 
if at the same time any effects associated with fmiteness of phonon ener- 
gies are also isotropic. Since it is known that the latter effects are impor- 
tant for incipient quantization effects in longitudinal magnetoresis- 
tance, 8 this last assumption deserves further investigation; however, we 
shall make it and proceed. Then, for the conductivity along H of all the 
valleys i in parallel, we have 

».<H) = E,[X/'V<(ff ( °) + X, (t) <r,(2? (f) )], (3G) 

where the coefficients \i %) , X t (l) depend only on the orientation of H 
with respect to the axis of valley i. Therefore, since ai = a t = do at 
H = 0, 

<r,,(H) 



<r(0) 



>i(fl (0 ) x, (0 , fr(fl (<) ) A/ n <:17) 

*» E(Xi w + x. w ) *o E(x/° + x t (<) ) 

i i 

Now, in the absence of quantization, a t is of order H~ or less as II — » °o } 
while a i remains finite; when quantization is just starting to be impor- 
tant and ull/c » ], a t will still be «* i . Under our assumptions, the 
dependence of both these quantities on H is given by the theory for the 
simple model with isotropic effective mass. If we keep only the first term 
of (37), we have the result italicized above. For some orientations, how- 
ever, the coefficient involving X's in the second term of (37) is larger than 
that in the first term by a factor of the order of a fraction of m\\*/m ± *. 
This explains the qualification in the statement of the theorem. 

3. In tlie approximation of isotropic scattering over an energy shell (i.e., 
n = Ti), any contribution to the transverse magnetoresistance ratio 
p(H)/p(0) which goes as H~ must for a cubic crystal be independent of 
the orientations of E and H, in the range y.H/c » 1, and any contribution 
which goes as H has only the modest orientation dependence shown in Table 
XIII. 

We interpret our assumption in the same way as for item 2 above, 
and with the same degree of risk. The proof is again quite simple if we 
use the transformations (34) and (35). Since E is covariant, H contra- 
variant, EH = implies EH = in every valley. When nH/c ^> 1, 
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Table XIII 
Orientation dependence of a hypothetical orbital-quantization term in 
the transverse magnetoresistance which is proportional to H. The value 
r = m* x /m\ = 0.052 has been assumed. 



nirection of 


Value of (40) = relative value of 


H 


j 




[001] 
[Oil] 

roil] 


J.H 

[100] 
[011] 


1.03 
1.74 
1.28 



E ~ RHj, and so, with j along x and H along z, 

pUH) = j- 2 a vy E°- = B*H\ y = R 2 H 2 £ *jW^— 

i & 



(38) 



where, as before, crt(H) is the conductivity of the simple-model spherical - 
ized valley transverse to a magnetic field H, and the summation is over 
the different valleys. In the absence of orbital quantization, at would go 
as H~ 2 , and so (38) would be independent of the magnitude of H. A 
term in p xx going as H can arise only from a term in & t independent of 
H. If bd t is such a term, the corresponding increment in p xx is 



5 PxI (H) = R-H'batY, 



E^_ 
E- 



(39) 



But the summation is a quadratic function of the components of E with 
the symmetry of the crystal, hence it is isotropic for any cubic case. This 
proves the first statement. A term in p xx going as H must come from a 
term in a t going as Jt~ , so the relative magnitudes of this term in various 
orientations must be proportional to the values of the quantity 



?(&)(?)'■ 



(40) 



These values are shown in Table XIII. 

4. Theory and experiment agree that the effect of orbital quantization on 
longitudinal magnetoresistance in 'pure germanium is very small as long as 
Tucr/kT < 1 or so. 

Calculations 8 based on the simple-model theory of Argyres 7 give p(H) 
actually about 15 per cent below p(0) for fuo c /kT = 1, rising to p(0) 
for hco c /kT ~ 2 and increasing linearly at high fields. The initial depres- 
sion, however, is due to the fact that, in this theory, the scattering is 
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describable by a relaxation time r(e) which goes to zero when e = 
(n + \)hu c . Since r _1 is proportional to the density of states in energy, 
the mean value of t _1 in an energy range several times fua e in width is 
close to the value calculated without quantization. Because of the 
fluctuations, the mean value of r over such a range is greater than the 
reciprocal of the mean of t _1 ; hence the conductivity is greater than in 
the absence of quantization, the current being carried preferentially in 
the energy ranges with the highest t's. Now these fluctuations in t will 
be considerably smoothed out if the energies of the phonons responsible 
for the scattering are even a fraction of fiu c , and crude numerical calcu- 
lations 8 show that a reasonable amount of such smoothing will change 
the maximum reduction of p from 15 per cent to a few per cent. Thus, 
we expect theoretically that the effect of quantization will be small for 
at least the range fto) c /kT < 1. 

The highest values of ho) c /kT which have so far been obtained experi- 
mentally for material with essentially acoustic scattering have been those 
of Furth and Waniek 12 for germanium in pulsed fields at room tempera- 
ture. Their maximum field of 460,000 gauss corresponds to hco c /kT = 2.1. 
Up to this field, their p(H) shows nice saturation, agreeing fairly well 
with the predictions of the theory ignoring quantization. By contrast, 
a combination of the latter curve with the quantization correction com- 
puted for zero phonon energies gives a curve which is rising conspicuously 
at the highest fields, and this is in marked disagreement with the obser- 
vations. Moreover, the saturation values of longitudinal magnetoresist- 
ance in both the [100] and the [111] directions agree with the predictions 
of the unquantized theory with w = m^T^m^Tx equal to about 0.057, 
a fairly reasonable value. 

5. In the region nH/c » 1, ftu c /kT < 1 and, for acoustic scattering, the 
fractional changes in Q P (H.) and in p(H) attributable to orbital quantiza- 
tion either should both be «./, or else should be roughly equal. 

Since the contribution of a phonon mode q to the Peltier flux is propor- 
tional to [c(q)] 2 T p /,(q), where c(q) is its velocity and r pA (q) is its relaxa- 
tion time, the statement just made will be valid if we can show that the 
crystal momentum lost by a current to the phonons is distributed in 
nearly the same manner among the different modes, regardless of whether 
quantization is assumed or ignored. 

It is easy to see qualitatively that, as huJkT becomes »1, the scatter- 
ing of the electrons must involve phonons whose wave-vector com- 
ponents _L_H become very large. This will undoubtedly cause the mean 
relaxation time of the phonons to become much shorter. Also, the change 
in the average direction of the wave vectors of the contributing phonons 
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will alter their mean relaxation time if the individual r p /i(q) are aniso- 
tropic. Finally, if quantization changes the relative contributions of 
transverse and longitudinal branches, the average [c(q)]"V p h(q) will be 
changed. We shall try to show, however, that all these effects are rather 
small when huJkT < 1. 

Since the theory of quantum magnetoresistance is rather complicated 
when TiUc/JcT < 1, we shall base our estimates on calculations of the 
distribution of crystal momentum among the different phonon modes for 
the case hu c /kT » 1. These calculations give the asymptote of the curve 
of (cTph)x\ against Tiio c /kT, an asymptote which should be very nearly 
valid for Tio) c /kT tt 3. Since the curve almost certainly is monotonic and 
comes in to its low-field limit with a horizontal tangent, one can esti- 
mate the order of magnitude of its deviation from the low-field value for 
hu c /kT < 1, or at least a rough upper bound to this. It will suffice for 
our purposes to give the argument for a semiconductor with isotropic 
effective mass, since we have shown above that, when the scattering is 
nearly isotropic over an energy shell, the transport problem for an aniso- 
tropic valley can be reduced to that for the isotropic case. 

For E parallel to H the appropriate quantum transport theory is well 
understood. A brief inspection of the formulas given by Argyres shows 
that, for TiuJkT » 1, the crystal momentum in the direction of the 
current which is delivered to the phonon modes of various wave numbers 
q is proportional to 



qieX »K8^kT }eXP 



-«y + fe a > i (4i) 



where we have taken the z direction to be that of E and H. From this, 
one can derive any kind of average which may be desired. Although our 
main interest is in an average t p a , hence in an average of something like 
q~ 2 or (f 1 with the weight (41 ), it will suffice to compute the more easily 
evaluated averages of qj, q v 2 , q z 2 , noting that the average of q~ will 
differ by a somewhat smaller percentage from its zero-field value than 
will that of q. It is clear from (41) that the averages of qj and q y 2 will 
go proportionally to w c , while that of q z ~ will be independent of a> c . The 
most convenient quantities to tabulate are the ratios of these averages 
to the corresponding averages for H = 0; such ratios are given for 
hcoc/kT = 3 in the first row of Table XIV. They differ by so little from 
unity that it seems safe to assume that, for ftta e /kT < 1, the average 
ct v h is very close indeed to the value at H = 0. 

It is worth remarking that the decrease of longitudinal resistance be- 
low the zero-field value, which the theory predicts 7,8 in the region 
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Table XIV 
Ratios of the average q x 2 , q y ~, q z 2 of the phonon modes, weighted in pro- 
portion to the crystal momentum they receive in the current direction, 
to the corresponding averages in the absence of a magnetic field. The 
ratios were evaluated from the high-field asymptotes for the value 
fioic/kT = 3. For the transverse case the predictions of two rival theories 
are given. Acoustic scattering has been assumed. 



Direction 


Theory 


(?i ! )Av. at S 
(?x 2 )Av. at o 

(linear in 11) 


(«1/ 2 )Av. at 3 

fol/*)Av. at o 
(linear in H) 


(?* 2 )Av. at J 
(<72 2 )Av. at o 


E 


H 


(constant) 


Z 
X 
X 


z 
z 
z 


Argyres 7 

Argyres 10 

Klinger-Voronyuk 9 


1.12 
1.12 
1.77 


1.12 
1.12 
3.37 


1.29 
-+0 



Tiu) c /kT < 2, probably has no very noticeable counterpart in the be- 
havior of the ((7x 2 )av etc. This decrease arises from the current being car- 
ried preferentially in the regions of energy just below the values (n + 
\)fwi e , where the relaxation time is longer than average. But the changes 
in the average q x 2 etc., for transitions at energy e as e passes through 
(n + \)fu»c are much less pronounced than those in the relaxation time. 

Similar calculations can be made from the theory of transverse quan- 
tum magnetoresistance, except that the correct form of the theory is 
still in dispute. 9,10 The alternative versions agree in predicting that for 
futiJkT » 1 the current should be carried predominantly in states of low 
h z , hence should involve transitions of low q z . If the levels had zero 
natural width and the phonons had zero energy, the integral for the 
current in the E-direction would diverge. In the last two rows of Table 
XIV, where we have tabulated the predictions of the two types of theory, 
we have therefore used the entry "— > 0" in the last column. The dis- 
tribution of crystal momentum in the current direction among the 
various modes turns out to be proportional, as far as q x and q u are con- 
cerned, to the last factor of (41) in the Argyres theory, 10 and to q* 
times this factor in the theory of Klingerand Voronyuk. 9 While the latter 
theory spreads the ^-distribution rapidly as H increases, the over-all 
average c~T p h reflects the combined behavior of all three components of 
q, and it seems likely that, for fuaJhT < 1, the average CT ph will not 
deviate from its zero-field value by more than a small fraction of the 
latter. 

By contrast, the values which these theories predict for the resistance 
p(H) are at least several times p(0) for fia> c /kT = 3, for a substance like 
germanium. 
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G. The high-field transverse magnetoresistance is so sensitive to inhomo- 
geneities, surface conduction and other perturbing influences that it cannot 
serve as a reliable indicator of the size of the orbital-quantization effect. 

At least, before it can be so used, many auxiliary investigations will 
have to be made to ensure the elimination of such perturbing influences. 
For the present purpose it will suffice to consider the case where the 
normal bulk conduction is in parallel with conduction by some other 
mechanism (e.g., surface conduction) which has a widely different Hall 
effect, say zero. Slight inhomogeneities in the bulk conduction itself can 
be shown to have a similar effect, 4 though they are probably not of 
major importance for our specimens. 

Let the .r-axis be along the length of the specimen, and let the mag- 
netic field H be in the z-direction. Let pb and p s be the effective resistivi- 
ties due respectively to normal bulk conduction alone and to the per- 
turbing conduction alone, and let Rb , R 3 = be the corresponding Hall 
constants. Then we have, assuming for simplicity that pb and p« are 
isotropic, 

jsy = -jby , (42) 

E v = Pb]/,,, — HbJJjbx = —Psjby , (43) 

E x = p,,jbx + RbHjb V = p„j sx . (44) 

From (43) and (44) we get the effective resistivity 

E x pM.Pb + Ps) + (R b H) 2 } 



PM = 



j hl + ./»*) (p h + Ps) 2 + (RbH)- 

* Pi 1 T 

L PbPs J 



(45) 



if p 3 » Ri,H » pi, . The physical meaning of (45) is that the Hall field 
due to the ^-conduction is partially shorted by the s-conduction, with 
consequent increase of the total dissipation. If p, is, say, 200 pb but 
RbH/pb is 7, the second term in the bracket in (45) will be 0.25, and the 
effect on a parabolic extrapolation to infinite field, such as we have 
made in Fig. (i, will be between two and three times this. 

The etching experiments described in Section V have shown that sur- 
face treatment can alter the observed transverse magnetoresistance by a 
very sizable amount at liquid-air temperature, expecially for H || [001], 
j || [100]. The subject merits a much more thorough investigation. How- 
ever, it should be emphasized that the high-field AQ effect is not com- 
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6 6 10 12 14 

MAGNETIC FIELD IN KILOGAUSS 



Fig. 17 — Values of the ratio of the transverse AQ(H) to Qp(0) for Sample 606 
at 20 5° and 44°K, showing the presumed effect of orbital quantization. The values 
of IQ P (0)| are 8600,xv/deg at 20.5° and 3730 at 44°. 

parably sensitive, either experimentally or theoretically; this is, of course, 
because there is no large Hall field in the AQ experiment. 

7. At 20° and U°K the effect of orbital quantization on the transverse 
Q p can be identified in the data, and it is small. Since the mechanisms of 
scattering of phonons are quite different at these two temperatures, these 
small values can hardly be due to an accidental cancellation of effects of 
quantization on p(H) and on CT ph , and it is probably safe to conclude that 
at liquid-air temperature quantization modifies p and Q p by only a few per 
cent. 

Fig. 17 shows the AQ data for Sample 606 at these two temperatures. 
Both curves start to saturate as H increases, then turn upward again. 
These upturns are undoubtedly due to orbital quantization; one may also 
suspect the presence of a quantization term varying less rapidly with H 
than the conspicuous part of the upturn, say linearly. One would like to 
separate AQ into a nonquantization part that saturates when fiH/c » 1 
and a quantization correction dependent on fiu c /kT. If, as the data sug- 
gest, a curve of the quantization correction alone against H would be 
concave upward, then we must surely underestimate the infinite-field 
limit without quantization if we assume it to equal the intercept on the 
H = axis of the tangent to the Q(H) curve at its most nearly hori- 
zontal portion. Thus, the difference between this intercept and the or- 
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dinate of the curve at 18,000 gauss ought to be an overestimate of the 
quantization correction. In this way we estimate that at 20.5°K, the 
orbital quantization correction is <10 per cent of Q P (0), and that at 
44° it is <3 per cent. The reasonableness of our partition of AQ into its 
two parts is attested by the values obtained for what AQ/Q P would be 
in the absence of quantization. Our intercept at 44° corresponds (if 
| AQ. | « 43/xv/deg) to a &Q P /Q P of 0.21, only slightly larger than the 
values found in the liquid-air range. The intercept at 20.5° corresponds 
to a AQp/Qp of about 0.33; this exceeds the high-temperature value by 
about the right amount to be attributed to the alteration of the phonon- 
phonon scattering law (see next paragraph). 

Now, at 20.5°K, the scattering of the phonons involved in Q p is almost 
entirely due to the boundaries of the specimen, while at 44° the phonon- 
phonon scattering exceeds the boundary scattering. Since these two 
types of scattering have quite different dependences on frequency and 
direction, it would not be reasonable to postulate that the effect of 
quantization on the average relaxation time of the participating phonons 
at both temperatures happens to be such as almost to cancel the effect 
of quantization on the magnetoresistance. The scattering of the elec- 
trons, incidentally, can be taken as approximately acoustic at both tem- 
peratures, since at 20.5° the observed Hall mobility, 215,000 cm 2 /vs, is 
about 0.8 of the estimated ideal value. We therefore conclude that, with 
acoustic scattering, the effect of quantization on the transverse p(H)/ 
p(0) is probably no more than 10 per cent or possibly 20 per cent at 
20.5°K. This conclusion is in violent contrast to the result of direct 
magnetoresistance measurements at this temperature, which for this 
sample (unetched!) gave p(H)/p(0) = 4.8 at 18,000 gauss. We believe 
the discrepancy is to be attributed to the great sensitivity of high-field 
transverse magnetoresistance to perturbing influences, especially surface 
conduction, as discussed in the Section V of the text. 

From the figures given above and the plausible assumption that 
above 44°K the orbital-quantization term decreases at least as fast as 
T~\ we conclude that, at liquid-air temperature and 18,000 gauss, this 
term is probably no more than a few per cent of the resistivity p(H), and 
almost certainly is no more than this for the thermoelectric power Q(H) 
or for Q P (H). It should be remembered, of course, that the effect of the 
quantization term on the extrapolation of p or Q to infinite field, made 
on a plot such as those of Fig. 6, will be proportionally larger than the 
effect at the field at which the extrapolation commences; for our para- 
bolic extrapolation procedure, the factor is a little under 2 for a term 
proportional to // and a little over 2 for a term proportional to H . 
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APPENDIX B. DERIVATION OF FORMULAS 

Iii this appendix we shall give some of the intermediate steps involved 
in applying the procedure outlined in Section III for the calculation of B 
and AQ at low and high magnetic fields. Our objective in doing so is 
partly to clarify the procedure used, but primarily to list a number of 
intermediate formulas which are especially useful or illuminating, and to 
give the final formulas in a form sufficiently general to be applied to 
other cubic many-valley structures (e.g., silicon), and to cases where 
the anisotropics are energy-dependent. 

Our starting point is the electron-group picture, with each ellipsoidal 
energy shell constituting a group g. The quantities we wish to calculate 
are the components of the Peltier tensor n, or, rather, the terms of order 
1, H, if in n at low H, and those of order 1, ir 1 at high H. For both 
cases we shall base our calculations on (6) or (9), which relate n to 
the partial conductivity tensors 6 g of the groups. These d , functions of 
H, are derivable from the transport equation (10) for group g, namely, 



-l 



■m*-d) ± 0dj X H - e 2 U^J E dn, (46) 

where the upper and lower signs, it will be remembered, are for electrons 
and holes, respectively. 

b.i Small H, Basic Formulas 

The conductivity tensor dd( = d„) of an energy shell is the coefficient 
of E in the solution of (46) for dy For H — > the usual way to expand 
dd in powers of H is by iteration. In a coordinate system oriented along 
the principal axes of the valley we find for the parts containing H°, H 1 , 
H 2 respectively: 

da J" = i (A) fe) dn, d*+ m = if * a ; (47) 

^-^(Sfe)^)?^* 1 ' (48) 

where 8 a p y is ±1 when a, 0, y are even (odd) permutations of 123, and 
is otherwise; and 

*„« = - (i 4 ) (A Z ( T? * ) Wp-Mi • (49) 
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b.2 Value o/QalH = 

Let us start by computing Q for H = in terms of the principal com- 
ponents II||(e), Tl ± (e) of n„ . From (6), (1) and (47) we have, for a 
cubic crystal with axially symmetrical valleys, 



1 q 2-^ \ "0<*t*' 

TQ(0) = ^Z a ii aa = « ■ 



J (0) 

ltd,,,. 



A* (0) 



I?/ 

<cn,,ri|/m,|*) + 2 (en^rx/m^*) 



(50) 



<6T ll An ll *) + 2Wm i *> ' 

where, as always, the angular brackets denote Maxwellian averages. This 
equation, of course, also holds for Q„ and Q p separately, with H||.a. — » 
n c or IIpii.j. . When the anisotropics p = H Pi /Ilp ± and w = m x *T\\/ 
ni\\*T X are energy -independent, the expression for Qp reduces to 

1^(0) - <£^ %+2"> . (51) 

<€Tj.) (2 + M>) 

When n p and -c depend on energy according to simple power laws, the 
Maxwellian averages occurring in these and subsequent expressions can 
be evaluated from Table XV at the end of this appendix. 

b.3 The Low-Field Value of B 

From (1), (5) and (6) we have, if H is in the z direction 

Tt(0) (1) 

TBH = lW» = [E iW 1 ^ P W + ?-$-. (52) 

If the crystal has over-all cubic symmetry, the evaluation of the first 
term on the right of (52) can be greatly simplified by noting that it is 
isotropic and may therefore be replaced by its average over all permuta- 
tions of the directions of x, y and z (keeping the system right-handed, of 
course). If this averaging is done before the summing on valleys, the 
result will be the same for all valleys. For a single valley it can be eval- 
uated by taking the coordinate axes along the principal axes of the 
valley. We find in this way, using (48), 

rr ivV"] = =f V*g*\<£±%> + <3^U> + Wilk)], m 

\_o Jxv 3c(e) _ m ± *- m\\*m ± * m\\*m ± * J 

where n is the total carrier density and, as usual, the angular brackets 
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represent Maxwellian averages. In the last term of (52) we have, from 

(ID, 

(1) phH 

» 

€ 



Pxy 



,(0) 



= ± 



(54) 



where p H is the low-field Hall mobility. We may use these and the ex- 
pressions of Ref. 3 — also obtainable from (47) and (48) — 



(0) 



= p 



co- 



ne 

m 



P-h = e 






+ 



2(«r 



2<er x ) 



+ 






m[|*m 



.1. 

rx>1 f2( e r x ) 
n±*J L ^i* 



+ 



mo* J 



(55) 



(56) 



to obtain finally 

Wiix) 



TB = =F 



Wj 



+ 



<«■! 



j.n x ) , <ct||t a ii 



W|| *m j. * 



I±nj!>l 
*ra x * J 



2<er 



Lm x * 2 rau*m x J 



(57) 



C 



This again applies to B p and J3 e separately, as well as to their sum, if 
IT||,j. and Q are given the appropriate subscript. When the anisotropics 
are independent of energy the formulas of Table I are easily derived from 
this with the use of (51). Note that the sign of B is independent of the 
sign of the carriers, since for electrons the upper sign is to be used in 
(57) with negative ITs and Q, while the reverse obtains for holes. 

All the formulas we have written thus far apply to any cubic material 
with axially symmetrical valleys. 

b.4 General Formulas for AQ at Small II 

For the low-field AQ we must take the second-order part of (9) and 
use (48) and (49). Alternatively, we can eliminate the dd° 's in favor of 
B. For the second-order part of the thermoelectric power in the a direc- 
tion we find, for extrinsic material, 



/„ (0) \ 



].. + Q 



(0) 



L P (o) _ 



-)} 



(58) 



*£ TB [Z* )\(H 2 -H a 



with the usual convention of upper sign for n type, lower sign for -p. 
The quantity p aa /p is just the magnetoresistance. Our task is now to 
compute the first term of (58) for various orientations. 
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Since in cubic material three phenomenological constants suffice to 
describe Q'"' for all possible orientations of H and YT, only three evalua- 
tions of (58) need be made. We can show further than when the valleys 
have [100] or [1 1 1] orientations in the Brillouin zone, their axial symmetry 
alone gives a further relation between the three constants, at least to 
the accuracy of the electron-group approximation, so that only two ori- 
entations need be computed. The argument, which is a generalization of 
the corresponding one for magnetoresistance, was given in Appendix C 
of our first paper. 1 In terms of the qb , q,- , q<i defined under Table IV the 
results are 

q, = -q b for [111] valleys, (59) 

q b + q* + 8i = for [100] valleys. (60) 

The evaluation of (58) can take a variety of forms, depending on 
whether p <2) and B are expressed explicitly in terms of the x's and n's, 
or are left as themselves. These forms can be converted into one another 
by use of (57) etc., and of 



(2) Z-, Vgaa 

Paa a 



z- 



p (0) o-IO) 

or 



(«I|H) 

(61) 



tE\ + 



IC'ffM 



(2) "I 



r (ll) 



(«±H). 



c 

The expressions for Zo °W <2> can be taken from the magnetoresistance 
literature, or (for [111] valleys) from (65) below or its longitudinal 
modification, by setting n = 1. 

b.5 AQ e at Small H with Energy-Independent Anisotropics 

As with the Nemst coefficient, the evaluation of Q e l2) for the case of 
an energy-independent anisotropy of x is a little easier than that of Q p . 
When H is parallel to the direction of measurement (the a direction), 
the last term of (58) drops out and, with n replaced by Te/e, the ratio 
of the first term of (58) to the second is a function only of the energy 
dependence of x. Using (61) we find easily, with p (2> = Ap and with the 
upper sign for electrons, the lower for holes, 

AQ e = Q,J 2) = =F? - r^l - £&1 -% (a longitudinal). (62) 
For transverse or general orientations of H we must use (57) to evaluate 
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the last term of (58) and must modify the relation of d„ z to p ; a 
little algebra gives 

AQ e V~a + 2 el\M(e) (ei*)(e)J p(0) 



These formulas are valid for any multivalley band structure in a cubic 
crystal. Still more general expressions for AQ e , valid for [111] valleys if 
the aiusotropy of t is not energy-independent, could be obtained from 
(G5) below. 

b.6 AQ and AQ P at S7nall H 

Because of the anisotropy of n p , the first term of (58) is a little more 
difficult to evaluate for the phonon-drag contribution or for the total 
AQ. We may again eliminate the sum over different valleys by noting that 
AQ is unchanged if we apply any operation D of the cubic symmetry 
group both to H and to VT, the a-direction in (58). If we average the 
contribution of each valley over all such D, the result will be the same 
for all valleys. If u (T) is the unit vector in the direction of VT, the result 
can be written, in the principal-axis system of some one valley, 

u -EDO,.* ].« --yi.^^v*'* - (64) 

where ( ) D means the average of the components shown over all sets, 
Z>H, Dn {T) . 

So far, the equations are valid for any cubic multivalley material. We 
shall now specialize to valleys along the [111] axes of the Brillouin zone. 
For this case and for H along a [100] axis, vT_|_H, (64) can be evaluated 
to 



,i< T > V Ttt a (2) 1 n (T) = - ng4g ' T? ( 6n ll r H 2T ^) _l I (eHimiTj. 2 ) 
u •Z.HVd, J-u - c2 [9 ( € ) m|| *2 mx * ^ 9 ( e ) Wl| * mx *2 



, 4 (cn x r||Tx 2 ) , 2 (eU ± rJ) l 
9 <€>m|*TO i « 9 (e)m ± * 3 J ' 



(65) 



When H and VT are both along a [100] axis we find an expression of 
similar form to (65) but with the coefficients f, £, £, f replaced respec- 
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tively by f , — ■§, — ■§, f. These expressions apply, of course, to either the 
electron or the phonon-drag contribution, or to their sum; one need only 
affix the proper subscript to all Il's. 

For the most general case the expression for AQ in terms of moments 
of II||,x is obtained by inserting (65), or its longitudinal modification, 
into (58) ; one may or may not wish to express the Q ° and B occurring 
there by (50) and (57). The magnetoresistance can be expressed if de- 
sired by (61) and (65) with n = 1, and p (0) and n„ by (55) and (56). 
The last entries in Table IX were obtained in this way. 

If we assume that p = n^/ripx and w = r\\m ±* / t ± m* are energy- 
independent, we get, after a little reduction, 

AQ T 1 = a fe*Yr, + (2 + «)• ((jJM _ <<r)W) V| (66) 

Q 'jm J '\ c ) V> ^ 3(1 + 2w)\ <«")« (MtW)}' (m " 

Jooi p \ C / 

2 (2 + w) 2 (l -V))(l - pw) («r) 2 (enr 3 ) 



= Q, 



3 (2 + pw)(l + 2u>y <enr)<€T 2 ) 2 

(vjJIY 2 (2 + to)(l -to) [ (l-w)(er)(eT 3 ) 
\ c J 3 (1 + 2wY L <«" 2 ) 2 

(2 + w)(l -pw) (er)\ ( 



(67) 



Hr 3 >] 



(2 + pw) <eIlT)<6f 

It is perhaps worth remarking that there exist a variety of checks 
which can be applied to calculations of this sort. For example, AQ should 
always vanish if all n„ are the same multiple of the unit tensor. Again, 
setting n a m*--e _1 should give the magnetoresistance. 1 These and other 
tests of the same sort are very helpful in uncovering algebraic errors. 

b.7 Large H, Basic Formulas 

The leading term in the expansion of the solution of (46) in powers of 
H~ l is especially simple when the geometry is such that E is either 
exactly parallel or exactly perpendicular to H. As these cases include all 
the orientations dealt with in this paper, we shall give first some formulas 
based on these specializations. 

When E and H are parallel, the part of dj normal to H must be of 
order EH~ X , since, if it were of lower order in H~\ the term of (46) in 
dj X H could not be compensated by anything else. Therefore, only the 
part dj B of dj parallel to H can be of order EH°. It is obtained by taking 
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the component of (46) along H: 

d J- - "' (n) iJ^X ■ ifE||H > (68) 

with neglect of terms of order H~ . 

When E and H are perpendicular, as they always are when H is per- 
pendicular to a symmetry plane and the total j is in this plane, we can 
start the solution of (40) by taking its cross product with the unit vec- 
tor u (/,) in the direction of H. The Hall field now dominates E, so all 
components of rfj are of order EH~\ and the component normal to H 
measures in essence the Hall velocity cEH~ x . We find, to the first order 
in H~\ 



(part of dj i H) = ± ec (j^\ u (//) X EH~ l dn if E _L H, 



(69) 



where the upper sign is for electrons, the lower for holes. The part of dj 
parallel to H is now obtained by taking the dot product of (46) with 
u (//) and using (69). The result of combining the two is, to the first 
order, 

4 - ± ec (£) [««- - (m ; n ;;, ) -; U "" ]| KB ± H, (70) 

where now we have used the label H X E in subscripts and superscripts 
to designate the corresponding coordinate direction, just as we used H 
above; u (HXE) is the unit vector in this direction. 

The expressions (68) and (70) will suffice for the calculation of the 
high-field limit of Q; to get the high-field behavior of the Nernst effect, 
however, for Table II, we need the conductivity tensor to one higher 
order. Rather than writing this down explicitly, we shall outline the 
result of a systematic solution of (46) by iteration of the procedures just 
used. Let us write for the conductivity tensor of the shell 

rfd(H) = dd <oo) + dd ( - i} + dd (_2) + • • • , (71) 

where the superscripts °o, —1, —2, etc. denote respectively terms of 
zeroth, first, second, etc. orders in H~\ With a corresponding notation 
for the components of dj, we find that, for fc > 1, the part of dj ( ~ k) 
perpendicular to H can be obtained in terms of dj (-fc+1) by taking the 
cross product of (46) with u (H) , and the part parallel to H from both this 
result and the dot product of (46) with u <fl) . The general term can be 
written concisely but abstractly in the form 
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d* ( - k) = - (jj?Xe 2 (£) (G-U-Sr-G-U-G^n, (72) 

if /< ^ 1; the direction of E is now unrestricted. Here S stands for the 
tensor m*-t -1 , while 

G a0 = 8 a0 - 8 n ~ l u m tm £ Wx <w) ^ , (73) 

x 

G + is the transposed tensor, and 

U* = EW"/"' (74) 

By use of the identity G + S = SG one can easily show that (72) is sym- 
metrical for even k, antisymmetrical for odd, as it should be. 

B.8 Limits of Q c and Q,, at Large H 

For either Q p or Q e , or their sum, we have, from (6) and the Kelvin 
relation, the high-field limit Q <00) given by 

TQ aa M = £ (iW V).. + Z Ov*.™ V +1) )« , (7.5) 

a ii 

where a labels the rectangular direction in which Q is measured, and 
where, by analogy with the notation already used, p <00) is the limiting 
resistivity tensor at infinite H and p (+1) is the high-field Hall tensor (12). 
When the a direction is parallel to H, the only non vanishing com- 
ponent of o-„ <00) is, by (68), the aa one. Moreover, in this case p0« (+1) 
vanishes for all /3, so only the first term of (75) remains. In this term we 
can replace p aa {x) by (£„ (r yaa <0O> ) -1 , so we have 

£ (enJ'Vdn*-^ 1 ) aa {i) ) 

TQj x> = -^ i — ts if H || a, (76) 

\-< , , I if. — 1\ (0\ ii J * ' 

2_w(m*-T ) aa ) 
i 

where i labels the different valleys. 

If the different components of t all have the same energy dependence, 
the relative magnitudes of the currents (68) for shells of different ener- 
gies will be the same as for H = 0, so, as (76) and (50) show explicitly, 

a (x) = Q.(0) parallel to H. (77) 

If the valleys have [111] orientations in a cubic crystal, the general 
form (76) becomes, for symmetry directions, 
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TQ (col 001 = (e[ll p| | + 2n p x][(m i ,Vr ll ) + 2(m x */rx)r 1 ) (?g) 

Wp Joox 3<e[(m l */r|) + 2(m i */Ti)r l > 

]on 
= 

, (79) 
(tIl pjL T ± /m x *) + (e[n px + 2n„||][(m A */Tj + 2(m B */rn)]" 1 ) 

<erxM x *) + 3(e[(wxVrx) + 2(m,*/ T| )]- 1 ) 

]in 

, (80) 

(en P ||Tn/m||*) + 3(e[n p |[ + 8n p J[(ffl ll Vr l |) + 8(CT A */rx)] ') 

(erii/m,,*) + 27<e[(m,*/r E ) + 8(m x *Ax)] _1 ) 

As usual, these equations also hold for Q e and Q p individually, if IT|| and 
n A are given the corresponding subscript. Further simplifications are 
possible if the anisotropics are energy-independent, as assumed in Table 
III, or if the t's are « e~~ , etc. 

When the a direction is perpendicular to H and H is normal to a sym- 
metry plane, (70) suffices for the evaluation of (75). For this case, the 
first term of (75) vanishes, since (r (00) has only an II H component and 
p (c0> has no aH component. We find easily 



TQ.™ = Ny- 1 ^ 1 Z ((nj fl - 1^) n„ ff (i) ) (« ± H), 



(81) 



where % runs over the different valleys, N v in number, and where, as 
above, S (l) is the tensor m*-x _1 for the tth valley. This holds for any 
multivalley structure when the symmetry is such that E is exactly _1_H. 
The entries of Tables III and IX are specializations of (81). 

B.9 Asymptotic Behavior of B as H — * co 

Let H be in the 2-direction, VT in the ^-direction. From (1), (5), (6) 
and the symmetries which we proved above from (72), we find 

r [Sl = 4? n »- ( -'""-" <+ " + - < " ,, - / " ) l/ (82) 

where, it will be remembered, the superscripts —2, —1, oo, +1 denote 
quantities going as If 2 , If 1 , H°, H as H — > oo. Here we may use (12) 
for p (+1) and get o- (_1) from (70). Since p +1) has only xy and yx com- 
ponents, the first term on the right of (82) involves just (^ fl n - o- £ , (-2)- ) „ . 
Use of the explicit evaluation of d„ (_2) = dd (_2) from (72) gives 
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(83) 



<*•*">-- -»£[(^ -¥■)(*• -¥0 

where, .as above, S is the tensor m** -1 for group g. If the directions of 
H and VT are such that p (oc) is diagonal in our coordinate system, the 
second term on the right of (82) involves (li g -(r ''~ 1) )xu > for which we 
find 

(n„.dd<-»)„ = ± c -|| l (n,„ - W.) , (84) 

where, as always, the upper sign is for electrons, the lower for holes. By 
the use of (50) and the relation 

(0) 

" = T p?-^ (85) 

ft R(co) 

we can evaluate (82) in the form, 
1 [ d(BH) ~\ (Mr 1 + 2S x -*))» 



V* / (lT CO Hxz' S zy \/r, («) S xi Syz \ 

/ (0 H i * S» * w g to ^f__\\ 

(eQS,-* + 2S ± ~ 1 )) v / / (l , _ TjJ»Sj»\\ 



(86) 



1 

where the superscript i labels the different valleys, N v in number. This 
equation applies to B e orB p or their sum, if all ITs, B and Q are given the 
appropriate subscript; it is valid for any multivalley structure with 
over-all cubic symmetry if z is normal to a symmetry plane and x is 
along a two-fold or four-fold axis in the plane. The entries in Table II 
result from tedious but straightforward specializations of (86) to [111] 
valleys and energy-independent anisotropics. 

b.io Maxwellian Averages 

As an appendix to this appendix we list here the Maxwellian averages 
of various powers of the energy, which are needed for the evaluation of 
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Table XV — Maxwellian Averages 
of Powers of Energy 



» 


«*)"> 


-1 


2 


3 

4 


1.38 


1 

3 


2tt- 1/2 


1 

4 


1.02 


1 


1 


I 

4 


1.04 


1 
2 


2w~ 112 


3 

4 


1.28 


1 


3 


!f 


1.81 


3 


iw-u* 


T 
4~ 


2.88 


2 


1 5 
4 


n 
5 


4.99 




127T- 1 ' 2 







the various expressions we have derived when the t's and IT's are as- 
sumed to have power-law dependences on energy. Since the number of 
Maxwellian carriers in the energy range e to e + de is proportional to 
e 5 exp ( — e/kT)de, we have for any quantity Z 



[ Ze h exp (-e/fcT) d* 

(Z) ^~ 

f e h exp (-e/kT) de 
Jo 



(87) 



The results for various powers of the energy are given in Table XV. 
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